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ABSTRACT 

When  investigating  the  detection  performance  of  a  passive 
homing  torpedo  used  against  shallow  draft  surface  ships, 
certain  environmental  factors  such  as  the  rough  sea  surface 
and  the  bubble-dominated  inhomogeneous  layer  near  the  sea 
surface  have  to  be  considered.   This  thesis  attempts  to  gain 
some  insight  into  the  behavior  of  a  homing  torpedo  system 
during  its  critical  attack  phase,  as  well  as  getting  some 
indications  of  the  relative  importance  of  the  scattering 
mechanisms  and  the  induced  tactical  limitations.   An  ideal- 
ized propagation  model  was  used  as  reference  of  comparison. 
For  a  given  sea  state  and  target  speed  the  results  stress 
the  importance  of  low  operating  frequency  as  well  as  a  high 
maximum  turn  rate.   They  also  point  to  the  importance  of 
having  a  search  depth  below  the  bubble-dominated  subsurface 
layer,  and  a  variable  speed  capability  during  the  torpedo's 
attack  phase. 
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I.   INTRODUCTION 

The  following  analysis  examines  several  factors  that 
limit  the  detection  performance  of  a  passive  homing  torpedo 
with  the  mission  objective  of  countering  shallow-draft 
targets  in  Norwegian  coastal  waters.   Generally,  these 
factors  can  be  divided  into  three  main  groups: 

-Environmental  factors  in  the  ocean 

-Electrical  mechanical  and  hydrodynamical  factors  in 
the  torpedo  system 

-The  users  tactical  situation 
The  factors  that  are  generated  in  the  ocean  itself  are  the 
principal  subject  of  this  analysis.   As  our  interest  is 
confined  to  the  layer  immediately  below  the  surface  of  the 
ocean,  later  called  the  subsurface  ocean  layer,  the  main 
factors  affecting  the  sound  propagation  are: 

-Scattering  and  absorption  due  to  the  bubble-dominated 
inhomogeneous  subsurface  ocean  layer. 

-Scattering  from  the  rough  sea  surface 
The  concentration  of  air  bubbles  and  the  roughness  of  the 
sea  surface  are  determined  by  the  windspeed.   The  effective- 
ness of  both  these  scattering  mechanisms  depends  on  the 
frequency  of  the  incident  wave  and  the  geometry  of  the 
source  and  receiver.   The  following  analysis  is  limited  to 
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high  frequencies  in  the  region  of  30-60  kHz  which  are 
characteristic  of  existing  torpedo  systems.   At  high  fre- 
quency and  low  grazing  angles  for  the  incident  and 
received  signals,  the  phenomena  of  "shadowing"  of  the 
surface  by  other  parts  of  the  boundary  occurs.   Under 
these  conditions,  the  effect  from  the  inhomogeneous  sub- 
surface ocean  layer  becomes  increasingly  important.   Ob- 
viously both  the  above  mentioned  scattering  mechanisms  are 
present  simultaneously.   Often  these  two  effects  cannot  be 
resolved  either  theoretically  or  experimentally,  as  any 
signal  with  a  finite  duration  will  be  scattered  from  the 
space  near  the  surface  simultaneously  with  that  from  the 
sea  surface  itself. 

In  order  to  adequately  describe  the  scattering  mech- 
anisms, this  analysis  starts  with  a  presentation  of  the 
oceanographic  background  material  for: 

-Typical  windspeed  and  wave  height 

-Typical  sound  speed  profiles 

-Density  and  distribution  of  air  bubbles  in  the  sub- 
surface layer 

-Statistical  description  of  the  sea  surface 
The  analysis  proceeds  by  separately  estimating  the  effect 
of: 

-Scattering  from  a  randomly  rough  surface 

-Scattering  and  absorption  caused  by  an  inhomogeneous 

subsurface  layer, 
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and  comparing  their  relative  importance.   The  method 
employed  for  these  estimations  is  an  approximation  that 
is  a  combination  of  both  ray  and  wave  theories.   Ray 
methods  are  used  to  follow  the  acoustic  signal  from  its 
source  to  the  vicinity  of  the  scatterer.   Wave  theory  is 
used  to  calculate  the  actual  scattering  process.   Finally, 
ray  theory  is  again  used  to  follow  the  scattered  signal 
to  the  receiver. 

An  idealized  propagation  model  consisting  of  an  iso- 
tropic stratified  medium  will  be  used  as  reference  of 
comparison.   This  model  is  founded  on: 

-A  noise  source  from  a  cavitating  propeller  blade. 

-The  operational  characteristics  for  a  square  law 
detector  (ROC-curves) . 

-A  transmission  loss  model  based  on  geometrical 
spreading  and  absorption  losses  in  homogeneous  sea 
water. 
In  this  analysis,  the  passive  sonar  equation  is  used  to 
predict  the  performance  of  the  homing  system.   The  detec- 
tion range  encountering  the  two  scattering  effects  will 
be  obtained  from  the  sonar  equation  and  compared  to  the 
detection  range  based  on  the  reference  model.   Thus,  the 
difference  in  ranges  at  which  the  homing  device  just 
acquires  the  target  with  and  without  scattering  is  a 
measure  of  efficiency. 
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II.   SCENARIO  AND  TORPEDO  RUN  GEOMETRY 

In  the  Norwegian  coastal  waters,  the  primary  mission  of 
a  torpedo  system  is  to  counter  an  amphibious  force  consist- 
ing of  escorts,  supply  ships,  and  shallow-draft  landing 
crafts.   Typical  characteristics  of  these  three  ship 
types  are  as  follows: 
Supply  ships: 

-Displacement  5000  tons 

-Length  100  m. 

-Draft  6  m. 

-Speed  15-20  kts. 
Escort  ships: 

-Displacement  2000-3000  tons 

-Length  85  m. 

-Draft  3  m. 

-Speed  35  kts. 
Landing  craft: 

-Displacement  1000  tons 

-Length  80  m. 

-Draft  2  m. 

-Speed  18  kts. 
In  order  to  simplify  this  analysis,  moderate  sea  states  (SS  3) 
are  assumed.   Since  "moderate"  wave  heights  of  two  meters 
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are  appreciable  when  compared  with  the  two-meters  draft 
of  the  landing  craft,  the  possibility  of  an  acoustic  torpedo 
impacting  the  target  at  a  depth  of  two  meters  is  very  remote 
without  the  use  of  an  influence  exploder. 

Two  relevant  search  and  attack  schemes  will  be  con- 
sidered.  These  are  illustrated  in  Fig.  1  together  with  the 
operation  of  the  influence  exploder.   From  target  valida- 
tion to  completion  of  terminal  attack,  the  torpedo  contin- 
uously tracks  in  the  azimuth  plane.   In  Case  A,  ascent  is 
inhibited  after  enable.   For  Case  B,  ascent  is  inhibited 
after  the  torpedo  reaches  terminal  attack  depth. 

An  assumed  attack  depth  of  six  meters  is  consistent 
with  the  activation  range  of  influence  exploders  and  is 
deep  enough  to  preclude  wave  or  "free  surface"  induced 
disturbances  of  the  torpedo.   Success  of  the  attack  depends 
primarily  on  the  availability  of  maintaining  azimuth-plane 
steering  to  within  a  short  horizontal  range  of  the  target, 
and  the  subsequent  operation  of  the  influence  exploder. 

Case  A  is  of  particular  interest  to  this  analysis,  as 
both  the  search  and  attack-depth  are  within  the  subsurface 
layer. 
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III.   OCEANOGRAPHIC  BACKGROUND  MATERIAL 

The  oceanographic  background  for  predicting  typical 
and  extreme  conditions  of 

-wind  speed 

-wave  heights 

-bubble  densities  and  distributions 

-ambient  noise  versus  self  noise 

-sound  speed  profiles 
are  outlined  in  detail  in  Appendix  A.   Even  though  most 
data  have  beneral  validity  for  Norwegian  coastal  waters, 
the  region  above  6  8""N  are  of  particular  interest.   Thus, 
a  typical  area  combining  open  as  well  as  confined  waters 
can  be  represented  by  "Andfjord"  at  70  N,  where  the  oceano- 
graphic conditions  can  be  related  to  the  weather  station 
"Andenes,"  see  Fig.  2. 

Figures  3  and  4  [Ref.  1]  show  average  windspeed  and 
the  occurrence  of  significant  wave  heights  as  a  function 
of  time  of  year  at  weather  station  "Andenes,"  respectively. 
The  bulk  of  data  is  centered  around  a  windspeed  of 
Beaufort:4-5  (11-21  kts)  and  SS:3-4  (significant  wave  heights 
1-2  m) .   Table  I  gives  the  relationship  between  SS ,  wind- 
speed  and  expected  significant  wave  heights. 
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TABLE  I 

RELATIONSHIP  BETWEEN  SEA  STATE  (SS) ,  WIND  SPEED, 

AND  EXPECTED  SIGNIFICANT  WAVE  HEIGHTS  (H  ) 

s 


Class 

Si 

gnif : 

Leant  Wave  Hei 

.ght 

Beaufort 

Number 

in 

m 

Sc 

ale 

SS 

Hs 

0 

0 

0 

1 

0  -  0.1 

1 

2 

0.1  -  0.5 

2 

3 

0.5  -  1.25 

3-4 

4 

1.25  -  2.5 

5 

5 

2.5  -  4.0 

6-7 

6 

4.0  -  6.0 

8 

7 

6.0  -  9.0 

9-10 

8 

9.0  -  14.0 

11 

9 

>14 

12 
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This,  together  with  the  low  probability  of  having  an 
amnibious  operation  occurring  in  high  sea  states  (SS>5) 
justifies  the  assumption  of  moderate  sea  state  with  wind- 
speed  in  the  region  of  12  kts  and  wave  heights  of  2  m. 

The  bubble  data  distribution  taken  from  Ref.  2  was 
obtained  in  the  area  "Tromso"  -  "Bjornoya"  during  the 
period  June-November  1978.   These  data  correlate  very  well 
with  a  larger  body  of  data  obtained  by  H.  Medwin  [Ref.  3]. 

Figures  5,  6  and  7  [Ref.  2]  show  the  density  of  resonant 
bubbles  as  a  function  of  depth  with  windspeed  as  parameter 
for  the   12,  38  and  120  kHz.   As  seen  from  these  figures, 
the  number  of  resonant  bubbles  are  an  increasing  function 
of  frequency  and  windspeed,  and  a  decreasing  function  of 
depth.   Below  a  depth  of  approximately  15  m,  the  number  of 
bubbles  is  negligible  for  the  windspeed  of  interest. 

The  effect  of  SS  (windforce)  on  the  ambient  noise  level 
is  given  in  Fig.  8  [Ref.  4].   Shallow  coastal  Norwegian 
waters  are  typically  5-10  dB  noisier  than  the  corresponding 
deep  water.   However,  great  variability  caused  by  local 
ship  traffic,  fishing  fleet  activity,  marine  life  and 
local  wind  conditions  makes  ambient  noise  prediction  diffi- 
cult in  these  areas.   This  means  that  for  accurate  modeling, 
ambient  noise  prediction  have  to  be  done  at  each  location 
as  its  level  is  both  site  and  time  dependent.   However, 
Fig.  8  shows  that  for  frequencies  higher  than  50  kHz,  the 
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effect  of  windforce  on  the  ambient  noise  level  decreases 
to  a  lower  bound  determined  by  the  thermal  agitation. 
Based  on  the  above  discussion  and  experience  related  to 
noise  levels  for  torpedo  systems,  the  self  noise  will  be 
assumed  to  be  dominant  through  this  analysis. 

Figure  9,  obtained  from  Ref.  5,  shows  that  the  sound 
speed  profiles  usually  encountered  in  the  area  of  interest 
results  in  extremely  difficult  sonar  conditions.   This  is 
illustrated  in  Figs.  10  and  11,  which  show  worst-case  ray 
path  derived  from  Fig.  9.   In  addition,  the  presence  of 
bubbles  in  the  subsurface  layer  causes  the  sound  speed  to 
be  a  function  of  frequency.   The  above  two  factors  may 
frequently  be  the  ones  limiting  the  detection  range  of  the 
torpedo.   These  effects  can  be  minimized  by  selecting  an 
appropriate  search  depth  for  a  particular  sound  speed 
profile.   In  addition,  for  Case  B  the  corresponding  curved 
homing  trajectories  in  the  pitch  plane  give  an  error  in 
apparent  range  to  target.   This  effect  will  normally  be 
taken  into  account  by  devising  appropriate  attack  logic 
which  is  outside  the  scope  of  the  present  analysis. 
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IV.   THE  PASSIVE  SONAR  EQUATION 

A  measure  of  efficiency  for  a  passive  homing  torpedo 
is  the  detection  range  obtained  by  solving  the  passive 
sonar  equation  for  broadband  noise: 

SL+lOlogB-TL (geom) -aR+DI-NL-DT=0  (1) 


where 


SL  =  spectral  level  of  the  broadband  noise 

radiated  by  the  target  (in  dB  re  luPa/lHz  at 
lm)  . 

R  =  detection  range  (in  m) . 

a  =  attenuation  coefficient  at  the  center  frequency 
(in  dB/m) . 

DI  =  receiving  sensitivity  (directivity  index) 
(in  dB  re  lyPa) . 

NL  =  noise  level  at  the  receiver  in  the  bandwidth  B 
(in  dB  re  lyPa) . 

DT  =  detection  threshold;  the  signal  to  noise  ratio 
at  the  transducer  output  required  for  a  de- 
tection probability  of  p   and  associated  false 


alarm  probability  p„,  (in  dB) . 
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V.   REFERENCE  MODEL 

A.  INTRODUCTION 

In  order  to  produce  the  reference  for  the  analysis 
the  sonar  equation  is  solved  assuming  ideal  free-field 
conditions,  a  simple  noise  source  model,  and  a  generalized 
square-law  detector. 

B.  IDEALIZER  TRANSMISSION  LOSS  MODEL 

Because  the  presence  of  refraction,  scattering,  and  of 
ocean  boundaries,  free-field  conditions  associated  with 
homogeneous  (isovelocity)  and  unbounded  medium  seldom 
exist  in  the  sea.   However,  as  a  basis  for  comparison, 
the  ubigutuous  spherical  spreading  law  plus  an  added  loss 
term  due  to  "normal  absorption"  can  be  used  as  a  reference 
model  for  measuring  the  effects  of  the  previous  mentioned 
scattering  and  absorption  mechanisms.   Thus,  the  reference 
transmission  loss  model  can  be  expressed  as: 

TL=201ogR  +  aR  (2) 

where  the  absorption  coefficient,  expressed  in  dB/m,  can 
be  obtained  from  Fig.  14  taken  from  Ref.  8. 
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C.   NOISE  SOURCE  MODEL 

1.   General  Characteristics  of  Noise  Sources 

Sound  is  generated  in  a  fluid  medium  by  any  process 
that  causes  an  unsteady  pressure  field.   Physically  processes 
that  can  cause  an  unsteady  pressure  field  include: 

-Pulsation  of  a  boundary  surface  of  the  medium 

-The  action  of  a  nonsteady  source  on  the  fluid 

-Turbulent  motion  in  the  fluid 

-Oscillatory  temperatures 
It  can  be  shown,  e.g.,  Ref.  10,  that  each  source  mechanism 
mathematically  corresponds  to  a  dominant  order  of  multipole. 
If  all  sources  are  of  such  a  nature  that  their  time  variation 
can  be  described  by  a  Fourier  Integral,  it  can  be  shown 
[Ref.  10]  that  the  Helmholtz  Equation  is 

72P^(x)  +  f£puU)  =  -4*fjx)  (3) 

3Q  (x)  52T  .  . 

+  V-F  (x)  -  -^± 


3t        gov  '    dx.dx. 

Term  1   Term  2    Term  3 
where  the  right  hand  side  describes  distributed  source  terms. 
The  terms  on  the  right  hand  side  of  Eq.  (3)  have  the 
following  interpretations: 

Term  1:  mass  injection 
Term  2 :  external  force 
Term  3:   turbulent  shear  stress 
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In  the  long  distance  and  long  wavelength  approximation, 
it  can  be  shown  that  the  mass  injection  term  gives  rise  to 
a  simple  source;  a  zero  order  pole  called  a  monopole.   The 
monopole  radiates  omnidirectional  with  no  angular  dependence 
At  large  distances  the  pressure  field  from  the  monopole 
radition  is  that  of  a  point  source.   Examples  of  this  are: 
-Pulsating  bubbles 
-Cavitation 
The  external  force,  in  the  long  distance  and  long  wave 
length  approximation,  is  associated  with  a  dominant  dipole 
which  has  a  cosine  directional  pattern.   Examples  of  this 
type  of  radiation  is  that  caused  by  the  vibratory  motion 
of  an  unbaffled  rigid  body. 

Radiation  from  turbulent  shear  stresses  is  charac- 
terized by  a  lowest  order  term  of  quadrupole  nature. 

The  efficiency  of  the  source  terms  decreases  with 
increasing  dependence  on  the  spatial  derivatives.   This 
can  be  understood  when  recognizing  that  wave  functions  of 
the  general  form  f (x-ct)  have  a  time  derivative 

l-^-f  (x-ct)  |  =  cft(x-ct)  (4) 

which  is  magnitude  c  (sound  speed)  greater  than  the  spatial 
derivative 

| ^-f (x-ct)  |  =  f  (x-ct)  (5) 

dX  X 
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Other  factors  being  equal,  the  radiation  from  an  external 

force  is  small  compared  to  that  from  mass  injection,  and 

that  from  turbulent  shear  stress  is  the  smallest;  therefore, 

moncpole  radiation  is  the  dominant  term. 

Propeller  cavitation,  when  it  occurs,  is  usually 

the  dominant  noise  source  for  any  marine  vessel.   Submarine 

and  torpedoes  often  operate  at  a  depth  great  enough  to  avoid 

cavitation.   Surface  ships,  on  the  other  hand,  generally 

have  well  developed  propeller  cavitation  with  the  result 

that  their  radiated  spectrum  from  5  Hz  to  100  kHz  is 

controlled  by  this  source. 

The  basic  phenomena  of  cavitation  combined  with 

propeller  hydrodynamics  give  the  fundamental  characteristics 

of  propeller  cavitation  noise.   An  excellent  qualitatively 

discussion  of  this  can  be  found  in  Ref.  9:Chs.  7  and  8, 

from  which  the  following  is  extracted: 

Propeller  blades  are  rotating  twisted  wings 
that  produce  hydrodynamic  forces.   Depending  on  operating 
conditions,  they  experience  cavitation  on  a  number  of 
different  places.   Of  these  there  are  three  prominent  types: 

-Tip  vortex  cavitation 

-Hub  vortex  cavitation 

-Surface  blade  cavitation 
In  addition  to  the  two  types  of  vortex  cavitation, 
there  normally  are  two  types  of  blade  surface  cavitation: 
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-Back:   driving  face 

-Front:   suction  surface 
Of  all  kinds  of  propeller  cavitation,  surface  blade  cavi- 
tation on  the  suction  surface  is  normally  the  most  noisy, 
while  hub  vortex  cavitation  is  the  least  noisy. 
2.   The  Noise  Source  Model 

Due  to  lack  of  recorded  and  available  noise  data 
from  the  target  in  question,  the  noise  source  has  been 
generalized  on  the  basis  of  the  following  discussion  and 
assumptions : 

The  noise  source  will  be  build  up  around  a  surface 
blade  cavitating  propeller  operating  in  a  good  to  poor 
wake;  surface  cavitation  will  be  assumed  to  be  dominant. 

D.  Ross  [Ref.  9]  has  developed  an  approximate 

theory  for  cavitation  noise,  where  dimensional  analysis 

is  combined  with  the  basic  results  of  cavitation  theory 

that  the  acoustic  pressure  is  proportional  to  the  product 

of  the  collapse  pressure  of  the  cavities  and  the  volume 

of  cavitation  produced  per  unit  time.   From  this  synthesis 

it  is  found  that  the  total  acoustic  intensity  varies  as: 

pbsD(U..)3   U.    U 
I  =  K,,  -2 **_[(  t  )(JL-  -i)*]  (6) 


where 


kti        2       U.  .   U. 
r        ti   ti 


r  =  distance  of  the  hydrophone  from  the  source 

K.  .  =  the  cavitation  inception  parameter 
ti 
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U  =  blade  tip  speed 

utj_  -   blade  tip  speed  for  inception  of  cavitation 

This  expression  shows  that  propeller  cavitation  noise  power 
is  proportional  to  the  total  number  of  blades,  b,  the  blade 
chord,  s,  and  to  the  propeller  diameter,  D,  and  is  a  function 
of  the  tip  speed  with  the  dependence  on  the  tip  speed  being 
the  strongest.   The  different  blade  surface  sections  where 
cavitation  exist  are  uncorrelated  and  the  radiated  noise 
is  treated  as  a  single  monopole  radiation  so  that  at  a 
distance  r>>a  (where  a  is  the  characteristic  dimension  on 
the  source  region)  the  radiation  is  similar  to  that  of  a 
point  source  with  no  angular  dependence. 

Submarines  and  torpedoes  with  centerline  propellers 
have  a  relatively  symmetric  inflow  condition.   Surface  ship 
propellers,  in  contrast,  operate  under  nonuniform  inflow 
conditions.   Circumferential  wake  variation  causes  the 
radiated  sound  to  be  amplitude  modulated  at  the  blade  rate 
frequency.   Furthermore,  slight  physical  difference 
between  the  blades  produces  modulation  at  the  shaft  rate 
frequency.   These  amplitude  variations  gives  a  very 
distinct  characteristic  to  the  radiated  noise  that  can  be 
used  for  classification  purposes  to  reduce  the  probability 
of  false  alarm. 

The  most  complete  source  of  data  on  surface  ship 
radiated  noise  are  measurements  made  during  WWII,  reported 
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in  a  compendium  issued  by  the  U.S.  Office  of  Scientific 
Research  and  Development  (OSRD)  in  1945  and  declassified 
in  1960.   When  these  data  are  examined  the  radiated  noise 
is  found  to  depend  on  tip  speed  and  the  number  of  propeller 
blades  with  little  dependence  on  the  other  variables.   For 
surface  ships  near  cruise  speed,  the  source  level  for 
frequency  over  100  Hz  can  be  written  as: 

SL=SL'+20-201ogf ;  f>100  Hz  (7) 

where 

f  =  frequency  in  Hz 

SL1 =  overall  source  level  in  dB  re  lyPa. 
The  overall  level  can  be  expressed  as: 

SL'=175+601ogU  /25  +101ogb/4  (8) 

where 

U  =  TrnD  (9) 

n  =  rotational  speed  (rpm) 

D  =  diameter  of  the  propeller  (m) . 
The  above  expressions  are  used  as  the  basis  for  the  noise 
model  with  the  following  input  data: 

n  =  300  rpm  for  maximum  cruise  speed  of  15  kts. 

n  =  180  rpm  for  a  cruise  speed  of  10  kts. 

D  =  2  m. 

b  =  5. 

The  resulting  noise  spectrum,  in  dB  re  lyPa  at  1  m, 
as  a  function  of  speed  in  kts.  are  tabulated  in  Table  II. 
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and  plotted  in  Fig.  15.   A  one  sigma  region  (5  dB  of  un- 
certainty) is  incorporated  in  Fig.  15. 

The  above  noise  spectrum  estimation  agrees  very 
well  in  the  high  frequency  limit,  with  more  recent  studies 
by  A.  Lovik  [Refs.  11  and  12].   Here  it  is  found  that  the 
cavitation  spectra,  both  theoretically  and  experimentally, 
can  be  divided  into  four  frequency  regions,  as  illustrated 
in  Fig.  16. 

Region  I  is  dominated  by  noise  at  the  blade  fre- 
quency and  its  harmonics.   The  emitted  sound  is  caused  by 
the  volume  variation  of  the  main  cavity. 

Region  II  starts  at  the  bubble  frequency,  which  is 
the  reciprocal  of  the  lifetime  of  the  main  cavity.   The 

mean  power  level  is  found  to  decrease  with  increasing  fre- 

,-2.5 
quency  as  f 

Region  III  is  an  intermediate  region. 

Region  IV  associated  with  the  shock  waves  starts 
at  the  mean  collapse  frequency  f  ,  given  by  the  mean 
collapse  time.   The  power  level  is  found  to  decrease  as 
f   ,  as  in  Eq.  5. 

The  number  of  gas  bubbles  in  the  water  have  a 
pronounced  effect  on  the  high  frequency  cavitation  noise 
from  the  propeller.   This  is  illustrated  in  Fig.  17, 
obtained  from  Ref.  12,  where  the  power  is  found  to  decrease 
as  much  as  4  0  dB  with  increasing  gas  content. 


30 


Scaling  laws  are  developed  [Ref.  12]  for  each  region 
based  on  a  series  of  models  and  full  scale  measurements. 
These  laws  depend  on  the  dynamic  pressure  induced  by  the 
propeller,  the  model  ratio,  and  the  gas  content  of  the  water. 
The  full  scale  measurements  were  performed  in  cooperation 
with  the  Royal  Norwegian  Navy  and  the  Marine  Institute  of 
Norway.   The  model  experiments  were  performed  in  the  largest 
cavitation  tunnel  at  the  Ship  Research  Institute  of  Norway. 

In  summary,  the  scaling  of  cavitation  noise  was 
demonstrated  to  be  a  useful  tool  in  predicting  a  full  scale 
cavitation  noise  as  shown  in  Fig.  18  [Ref.  12]  which  compares 
measured  noise  spectra  for  the  model  and  full  scale 
measurements . 

For  the  high  frequency  region,  the  source  levels 
are  of  the  same  magnitude  as  predicted  by  the  WWII  empirical 
formula. 

D.   PASSIVE  MODE  RECEIVER  CHARACTERISTICS 
1.   Assumptions 

For  receiver  characteristics  assume  a  square  law 
detector  with  a  center  frequency  f-60  or  30  kHz  and  a  band- 
width B.   The  detection  scheme  is  shown  in  Fig.  16.   The 
principal  assumptions  employed  in  the  derivation  are  as 
follows : 

-Gaussian  signals  in  Gaussian  noise 
-Frequency  independent  signal  and  noise  spectra 
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-Integration  time  T  is  sufficiently  long  to 
permit  application  of  the  central  limit  theorem. 
2 .   Derivations 

The  detector  input  r(t)  is  assumed  to  be  a  zero- 
mean  Gaussian  process  composed  of  noise  alone  or  signal 
plus  noise  expressed  by  the  two  well  known  hypotheses 

HQr(t)  =  n(t)  (10) 

f^rCt)  =  s(t)+n(t) 
where : 

n(t):   noise  signal 

s (t) :   signal. 
The  two  signals  s(t)  and  n(t)  are  assumed  to  be  independent 

Assume  that  the  spectral  shape  of  s(t)  and  n(t) 

are  the  same,  such  that  H   and  H,  only  differ  in  the  total 

o      1    J 

power  level.   Then  the  detector-smoother  have  the  form: 


rlt) 

V 

iOttt 

J 

x(t) 

T 
0 

K 

In 

Schematic  of  detector-smoother 


and 


x(t)  =  r2  (t)  (ID 

Furthermore,  let  the  noise  variance  be  normalized  to  unity 
(for  convenience)  and  the  signal  variance  by  denoted  by  a2 

Var[n(t)]=l  (12) 

Var[s(t) ]=a2 
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Because  of  the  assumed  similarity  in  the  spectral  shapes, 
the  autocorrelation  functions  are 

Rn(T)=F*1[N(f)  ]  =  p(t)  ;  N(f)  is  the  noise  power   (13) 

spectral  density. 

Rs(x)  =  p(x)a2 

Rn(x)=  Hq  p(x) 

H1  (l+a2)p(i) 
Furthermore,  assume  that  the  integration  time  T  is  long 
enough  so  the  central  limit  theorem  holds,  implying  that 
K  also  is  a  Gaussian  random  variable. 

This  yields  that  the  probability  density  function 

of  the  output  variable  and  hence  the  detection  and  false 

alarm  probabilities  are  completely  determined  once  the 

mean  and  the  variance  of  K  are  derived. 

If  a  process  V(t)  is  wide-sense  stationary,  then 

T  T 

E[V(t) ]=E^  A/(t)dt  =i  /E[V(t)]=v  (constant)    (14) 
0  0 

Thus,  assuming  that  r(t)  is  a  wide-sense  stationary  process. 

E[K]=E[x(t) ]=E[r(t)]=l+a2  (15) 

and  similarly 

TT 
Var[V]=E[V2]-{E[V]}2=   -2-      ff    E  [V(t)  V(s)  ]  dtds-v2     (16) 

00 


TT  TT 


Var[V]=-^2-    ff  [Rv(t-s)-v2]dtds=^2-    ff    Covv(t-s)dtds 
00  00 
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where  Covv(t-s)  is  a  covariance  function. 
Then  letting 

x  =  t-s 

C   =   t+s,    and   substituting      into   Eq.     (7)    yields 


2T 


d<;dx 


Var[V]=^/yCovv(T) 

o 


T       V 


rar[V)=i   /  [i-M]Covv(T)dT 


(17) 


-T 


Consequently 

T 
Var[K]=i  j"   [l-M]Covv(t)dT 


(18) 


-T 


Then,  evaluating  the  covariance  function  from  the  auto- 
correlation function 

R  (T)=E[x(t)x(t-T)]  =  E[r2(t)r2(t-T) ]  (19) 

Since  r(t)  is  Gaussian,  the  above  fourth  moment  can  be 
expressed  as  product  and  sums  of  second  moments: 


R  (i)=r2(t)  »r2(t+T)  +  2r(t)r(t-T)  -r(t)r(t-x) 

A 


=R2 (0)  +  2R2 (x) 
r        r 


2\  2 


R  (x)=(l+a^) ^+2(l+az) 

A 


02(T) 


(2o: 


Thus,  the  covariance  function  is 
Cov  (t)=R  (x)-(E[x(t)]}2 


=(l+a2) 2+2(l+a2) pz(x)-(l+a2) 


2\  2 
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Covx(x)=2(l+a2)  V(t)  (21) 

Inserting  Eq.  (12)  into  Eq.  (9)  yields 

T 
Var[K]=i  y*  {[l-Ill]  2(l+a2)2p2(T)}  dx 
-T 

Var[K]-2<y  /  [l-M]  P2<x)dx  (22) 

-T 
If  T  is  large  compared  to  the  correlation  time 

TB>>1 
then  we  can  substitute  the  limit  for  Eq.  (22)  by 


00 

2  \   2 


Var[K]=      2(J;+°    }         /     p2(T)di  (23) 


—  CO 

oo 

/  XT  2 


2(1+0')  <        /     NMf)df 
T  J 

—  00 

If  we  further  make  the  assumption  that  the  signal  and  noise 
have  ideal  flat  bandpass  spectra: 

N(f)=(l/2B,  f-B/2<f<f+b/2  (24) 

0,  otherwise, 
Inserting  Eq.  (24)  in  Eq.  (23)  yields 

Var[K]=  (1+^)Z  (25) 

The  probability  density  function  for  the  output  of  the 
detector  have  the  following  form 
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pMH0) 


pfwIH.) 


=ThreshoU 


H0  choosen  - — 


— *-  H,  cl-ioosen 


Here  w  is  the  outcome  of  all  possible  signals.   The  false 
alarm  probability  is  obtained  by  integrating  the  conditional 
probability  p(w|H  )  over  the  outcome  space  for  which  to 


choose  H- . 


'FA 


p(w|HQ)dw  =  Q(-^-J 


(26) 


Similarity,  the  miss  probability=l-detection  probability 

is  obtained  by  integrating  the  conditional  probability 

p(w|H, )  over  the  outcome  space  for  which  to  choose  H  . 
1  1  o 


PM  =  ^D  = 


p(w|H1)dw  =  Q( 


(27) 


Further  defining  the  input  and  output  signal-to-noise 

ratios  as 

a2     3 
S/N (input)  =  y"  =  a 


(28) 


S/N (output)  = 


(Ui-U0) 


(29) 


Ol 
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For  the  square  law  detector  where 

Po  =  E[n2(t)]  =1,   p,  =  E[r2(t)]  =E[K]  =l+a2   (30) 
a02  =  Var[n2(t)]  =  i,   cr^  =  var[r2(t)] 


2 


=  Var[K]  =  ii±2li 


The  output  signal-to-noise  ratio: 


S/N (output)  =  (1+°2-1)2  =  BT  a"  (31) 

rU-t-a2)2,    (1+a2)2 
L   BT   J 

The  input  signal-to-noise  ratio: 

S/N (input)  =  a2  (32) 

The  probability  of  false  alarm: 

PFA  =  Q[/BT (j-1)]  (33) 

The  probability  of  detection: 


P   =  i-q[/bt  (1+q2;^)] 

1+a 


(34) 


The  Equations  (32)  ,  (33)  ,  and  (34)  are  plotted  for  a 
variety  of  realistically  encountered  conditions.   Figure  17 
gives  BT  versus  S/N  (in)  (Identical  to  the  detection 
threshold  DT)  for  various  combinations  of  p_  and  p__ . 
Figure  18  gives  p  versus  S/N (in)  for  various 
combinations  of  the  threshold  j  and  PFA«   This  constitutes 
the  ROC-curves  for  the  generalized  square  law  detector. 
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E.   REFERENCE  DETECTION  PERFORMANCE 

1.  Introduction 

As  both  the  scattering  and  absorption  are  frequency- 
dependent,  it  is  necessary  to  have  reference  models  for  both 
60  and  30  kHz. 

2.  60  kHz  Case 

The  range  dependent  portion  of  the  sonar  Equation  (1) 
-201ogr-aR 

is  plotted  in  Fig.  19  for  a=0.021  dB/m  taken  from  Fig.  14. 

Figures  15,  17  and  19  are  then  used  to  estimate  the 
detection  range: 

(a)  Select  p  ,  PFA/  B,  and  the  integration  time  T. 
The  detection  threshold  DT=S/N (input)  is  then 
found  from  Fig.  17. 

(b)  Select  the  speed  of  the  target  and  find  the  SL 
from  the  noise  source  model  (Fig.  15) .  Then, 
reasonable  values  for  the  receiver  sensitivity 
DI  and  the  self  noise  level  NL  yields  the  left 
hand  side  of  the  sonar  equation  (1)  except  for 
range  dependent  term. 

(c)  Use  Fig.  19  to  solve  the  passive  equation 
for  R. 

A  realistic  example  may  illustrate  the  above  procedure. 
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(a)   Entering  Fig.  17  with 
PD  =  0.5 

Pfa  =  10~6 

B  =  4500  Hz 
T  =  100  msec 


yields 


DT  =  -6.5  dB. 

(b)  Selecting  a  target  speed  of  12  kts.  gives 
(from  Fig.  15)  SL=100  dB.   Selecting  a 
typical  transducer  sensitivity  DI  =  -180  dB. 
Assuming  the  NL  to  be  dominated  by  self  noise 
of  typical  value  NL=-124  dB.   This  yields 

SL+101ogB+DI-NL-DT=87  dB. 

(c)  Figure  19  with 

-201ogR-0.02R=-87 
yields 

R=1200  m.  for  a=0.02  dB/m  and  f-50  kHz. 
The  influence  of  different  design  parameters  like  self 
noise  and  detection  threshold  on  the  passive  detection 
performance  is  now  easily  investigated  by  the  above 
procedure . 

Although  outside  the  main  scope  of  this  analysis, 

the  above  statement  can  be  confirmed  with  an  example.   As 

—  6 
seen  from  Fig.  18,  a  probability  of  false  alarm  p   =  10 

implies  a  threshold  setting  =  0.9  dB  above  the  noise  level. 
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Due  to  the  variability  of  the  noise  level  together  with 
the  practical  difficulty  in  accurately  setting  the  threshold, 
a  more  realistic  goal  for  the  threshold  would  typically  be 
3  dB.   Going  into  Fig.  17  shows  that  the  corresponding  value 
for  DT  for  pD  =  0.5  is  DT  =  0  dB,  with  a  corresponding  low 
value  for  the  p__ . 

Letting  DT=0  dB  and  keeping  the  previous  assumed 
values  of  SL,  DI ,  and  NL  yields: 

SL+101ogB+DI-NL-DT=80.5  dB. 

The  corresponding  detection  range  is: 

R  =  1000  m,  for  a  =  0.02  dB/m  and  f  =  60  kHz. 
Thus,  this  change  in  threshold  setting  caused  a  decrease 
in  detection  range  from  1200  m  to  1000  m  in  return  of  a 
significant  decrease  in  the  false  alarm  probability. 
3.   30  kHz  Case 

In  order  to  estimate  the  reference  detection  range 
for  an  operating  frequency  of  30  kHz,  we  utilize  the  sonar 
equation  (1) . 

Assume  that  the  receiver  has  the  same  generalized 
passive  detector  characteristics  as  in  the  60  kHz  case: 

DI  =  -180  dB 

NL  =  -124  dB. 

—  6 
DT  =  -6.5  dB,  based  on  pD  =  0.5  and  pFA  =  10 

However,  the  empirical  equation  (7)  for  the  source  level: 

SL=SL'+20-201ogf 
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shows  that  the  source  level  falls  off  as  f"2.   If  the 
dynamical  and  dimensional  parameters  of  the  propeller  are 
the  same,  SL  will  increase  by  +6  dB  when  the  frequency  is 
reduced  from  60  to  30  kHz. 

A  source  level  of 

SL=100+5=106  dB 
gives  a  range  dependent  solution  of  the  sonar  equation 

-201ogR  -aR=-87-6  =  -93  dB. 

A  plot  of 

-201ogR-aR 
is  given  in  Fig.  20  for  an  absorption  coefficient  a=0.01  dB/m 
taken  from  Fig.  14. 

Figure  20  then  gives  a  detection  range  of 

R  =  2400  m. 
Thus,  as  seen  from  these  ideal  reference  calculations, 
halving  the  frequency  gives  a  higher  source  level  and  a 
lower  absorption  loss,  resulting  in  a  doubling  of  the 
detection  range. 


41 


VI.   THE  EFFECT  OF  SURFACE  SCATTERING 

A.   OCEANOGRAPHIC  DESCRIPTION  OF  THE  SEA  SURFACE 

The  roughness  of  the  sea  surface  is  the  essence  of 
the  scattering  mechanism.   Thus,  to  adequately  describe 
the  scattering  of  sound  from  a  randomly  rough  sea  surface, 
it  is  necessary  to  formulate  a  suitable  description  of  the 
sea  surface  from  an  acoustical  propagation  point  of  view. 
Generally  the  shape  of  the  rough  sea  surface  is  most 
appropriately  described  in  terms  of  time  and  spatial 
dependent  random  variables.   However,  observation  of  the 
ocean  under  the  same  environmental  (meteorological)  condi- 
tions indicates  that  the  roughness  is  the  same  over  large 
areas  and  for  periods  of  at  least  several  hours.   The 
random  processes  responsible  for  the  structure  of  the  sea 
surface,  therefore,  can  be  considered  stationary  at  least 
over  periods  of  hours.   With  this  assumption,  the  sea 
surface  can  be  described  in  terms  of  the  statistical 
description  of  the  surface  displacement  function,  and  the 
distribution  of  signals  reflected  from  the  sea  surface 
can  then  be  related  to  this  probability  distribution 
function. 

The  most  significant  statistical  parameters  describ- 
ing the  scattering  mechanism  from  the  randomly  rough 

surface  are: 
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-the  mean  square  slope 
-the  mean  square  surface  height 
-the  correlation  length. 
Optical  measurements  made  at  sea  by  C.  Cox  and  W.  Munk 
[Ref.  13]  showed  that  the  sea  surface  with  an  arbitrary 
wide  continuous  spectrum  of  waves  is  characterized  by  a 
Gaussian  distributed  surface  slope.   The  mean  square 
slope,  determined  from  these  optical  measurements  is 

<Cl2>  =  Z*    =  (3+5.12w)xl0~3  (35) 

where 

W=wind  speed  in  m/s  measured  41  ft  (12.5  m)  above 
the  sea  surface. 
The  Gaussian  distribution  of  surface  slopes  implies  that 
the  surface  displacement  function  can  be  described  by  a 
Gaussian  probability  density  function  with  zero  mean 

<C>  =  0 
and  variance 

<C>  =  a2 
and  Gaussian  correlation  function 

<s(t)e(t+T)> 

The  mean  square  height  a2,  is  obtained  by  integrating  the 
frequency  spectrum  of  the  fully  developed  sea.   The  fre- 
quency spectrum  G  is  given  by  the  commonly  accepted  semi- 
empirical  expression  of  Piersom-Moskowitz  [Ref.  14]  as 
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fi 


G(«)  =  ^_  exp[.e  _or  (36) 


where 


Q,  =  frequency  (in  s   ) 
a  =  8.1 
S  =  .74 

QQ=   g/W  (in  s"1)  . 

W = wind  speed  in  m/s  at  19.5m  above  the  sea  surface 
g = gravitational  acceleration  (in  m/s2). 
This  gives 

00 

a2  =  /  G(Q)an   =  ^—  (37) 

o  4  6g2 

For  a  Gaussian  autocorrelation  function  expressed  as 

iMt)  =  —  <;(t);(t+T)>=  e~jZ/ rZ  (38) 

a2 
where 

T  =  correlation  length. 

The  following  relationship  for  the  mean  square  height  holds 

for  sea  of  small  roughness 

Z2  =  22—   ,  see  later  Eq.  (60)  (39) 

T2 

and  the  correlation  length  is  thus: 

T  =  /2  |  (40) 
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B.   SCATTERING  THEORY 

All  real  boundaries  are  rough  for  radiation  with  short 
enough  wavelength,  and  the  apparent  roughness  depends  on 
the  "viewing"  conditions.   The  wave  reflected  by  a  plane 
surface  has  the  same  properties  as   the  incident  wave  since 
the  radiation  is  scattered  coherently  and  there  is  a  definite 
relation  between  the  incident  and  scattered  waves. 

A  randomly  rough  surface,  however,  such  as  the  wind 
generated  ocean  surface,  scatters  radiation  in  all  directions, 
i.e.,  an  illuminated  area  is  visible  from  any  direction. 

Heuristically  there  are  two  distinct  approaches  to  this 
phenomena . 

1.  If  the  boundary  is  rough  most  of  the  radiation  is 
scattered  and  there  is  little  transmission  in  the 
specular  direction.   Thus,  the  attenuation  caused 
by  the  irregularities  can  be  included  in  the  trans- 
mission equation. 

2.  If  the  surface  is  truly  smooth,  it  can  be  assumed 
that  the  effect  of  the  boundary  is  to  supplement  the 
original  pressure  field  by  an  out-of-phase  image 
contribution.   For  a  randomly  rough  surface  the 
reflected  sound  neither  completely  cancels  the 
direct  sound  nor  adds  to  give  +6  dB  pressure  peaks 
of  the  interference  pattern.   For  a  rough  surface, 

this  supplement  is  a  small  fraction  of  the  direct 

path. 
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The  second  approach  will  be  used  with  the  simplified 
assumption  that  the  sea  below  the  surface  has  an  isotropic 
statistical  description;  i.e.,  the  mean  acoustic  velocity 
and  the  mean  density  are  assumed  to  be  constant  and  have 
negligibly  small  mean  square  fluctuations. 

The  estimation  of  the  scattering  is  based  on  an  approxi- 
mation method  employing  both  ray  and  wave  theory.   Ray 
methods  are  used  to  follow  the  acoustic  signal  from  the 
noise  source  to  the  vicinity  of  the  sea  surface.   Then, 
wave  theory  is  used  to  calculate  the  scattering  process. 
Finally,  ray  theory  is  used  to  follow  the  scattered  signal 
to  the  receiver. 

The  geometry  is  given  in  Fig.  24a.   The  origin  of  the 
coordinate  system  is  at  the  center  of  the  illuminated  area. 
The  x-y  plane  coincides  with  the  mean  of  the  rough  surface 
as  averaged  over  the  illuminated  area. 

The  source  and  receiver  are  at  distances  R,  and  R2 , 
respectively,  from  the  origin.   R,  is  the  xz-plane  and 
makes  the  angle  9,  with  the  z-axis.   R2  makes  the  angle  92 
with  the  z-axis  and  the  projection  of  R2  on  the  xy-plane 
has  an  angle  9~  relative  to  the  x-axis. 

For  high  frequencies  R,  and  R~  are  much  larger  than 
the  acoustic  wavelength.   Then  both  the  incident  wave 
and  scattered  waves  can  be  treated  as  nearly  plane  waves. 
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The  formulation  of  the  scattering  problem  will  be 
based  on  the  Helmholtz  integral  which  requires  known 
values  of  the  normal  derivatives  of  the  incident  and 
reflected  waves  on  the  boundary.   These  are  estimated 
by  using  the  neuristic  Kirchhoff's  approximation,  which 
assumes  that  the  wave  is  locally  reflected  by  a  plane 
surface;  i.e.,  an  approximation  restricted  to  a  surface 
not  too  rough  and  not  shadowed. 

Further,  the  receiver  derictivity,  as  indicated  by 
Fig.  24b,  will  be  used  to  limit  the  surface  area  that  is 
illuminated. 

Since  this  procedure  is  based  on  a  detailed  develop- 
ment by  I.  Tolstoy  and  C.  S.  Clay  [Ref.  15],  only  the 
main  points  will  be  outlined  here  to  bring  out  the  assump- 
tions made  and  the  inherent  limitations  of  this  approach. 

The  development  starts  by  considering  the  inhomogeneous 
wave  equation  of  the  general  form: 

_       1  32p(xxt) 

V2p(Xjt)  - =  -4tt  f  (Xlt)  (41) 

c2    dt2 

where 

f(xtt)  =  is  a  known  source  distribution. 
The  development  is  based  on  the  following  initial  assumptions 

-The  medium  is  homogeneous. 

-The  medium  is  bounded  by  some  surface  S,  onto  which 
an  incident  wave  impinges. 

47 


-The  boundary  is  characterized  by  the  specific 

acoustic  admittance  and  the  shape  of  the  boundary 
-The  incident  wave  is  harmonic. 
The  assumed  harmonic  source  implies  that  f(x,t)  can  be 
decomposed  into  a  Fourier  integral.   Furthermore,  assuming 
that  the  solution  of  Eq.  (41)  can  be  decomposed  in  time, 
we  arrive  at  the  Helmholtz  equation 

V2p  (x)  +  ^-  p  (x)  =  -4TTf  (x)  (42) 

UJ  Q  JJ  CO 

It  should  be  noted  here  that  L.  Fortuin,  in  [Ref.  16] 
showed  that  the  Helmholtz  equation  is  not  exactly  correct 
for  a  medium  with  a  time  dependent  boundary.   The  equation 
can,  however,  be  used  with  a  good  approximation  when  the 
time  derivative  of  the  surface  elevation  is  much  smaller 
than  the  speed  of  the  waves  through  the  medium.   For  under- 
water sound  waves  scattered  by  the  rough  sea  surface,  this 
means  that  the  wind  speed  has  to  be  much  less  than  the 
sound  speed;  a  requirement  easily  fulfilled  for  our  inves- 
tigation. 

Green's  method  allows  the  solution  of  this  linear  in- 
homogeneous  wave  equation  to  be  expressed  in  the 
heuristic  Helmholtz  integral  form: 

r  !  /     3p  (xl) 

p    (x)=/f    (xMG    (xxMdV+^UtGUx1)  — (43) 

w        v    u  w  4Y  an1 

8G    (xx1) 

_  p    (x:)  — Ma1 

an1 
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p. 


{G(x2xM  -^  -  p  (xM^-jda1      (44) 
an1    w    an1 


The  first  integral  on  the  RHS  of  Eq.  (43)  contains 
the  sound  sources  and  the  bulk  (volume)  scattering.   The 
second  integral  represents  the  surface  scattering  and  is 
taken  over  all  finite  surfaces. 

Now,  disregarding  the  direct  path,  the  signal  at  the 
receiver  is  given  by  the  surface  integral  alone: 

As), 

l*2;=4T 

s 

where  the  subscript  s  denotes  the  scattered  field.   In 
order  to  solve  Eq.  (44)  the  following  must  be  done: 

-Give  an  approximate  expression  for  the  incident 
wave. 

-Find  an  appropriate  Greens  function. 

-Make  an  approximation  for  p  (x1)  and  9p  /3n  at 

00  00 

the  surface. 
As  we  already  have  assumed   a  simple  harmonic  source, 
the  incident  wave  can  be  expressed  as: 

p'^^ft^e^^fe1"  (45) 

roo  2tt    R         K 

where 

IT:   power  output. 

D:   illumination  function. 

Assuming  kR>>l,  i.e.,  that  the  distance  of  the  source 
is  large  compared  to  the  wavelength,  the  wave  in  the 
bounded  ensonified  area  can  be  considered  as  a  plane  wave 
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characterized  by  its  propagation  vector 

-»■ 
x. 

k.  =  -k-  x 


1        I"*"  I 

thus  obtaining  the  expression  for  the  incident  wave: 

P  (i)(;  ,  =  bd  eikRl  eiiti  -5' 

Further  assuming  that  the  receiver  is  sufficiently  far 
removed  from  the  scattering  area  with  the  rest  of  the  sea 
surface  not  contributing,  then  the  scattering  area  acts  as 
a  small  induced  source  in  a  free  space  and  we  can  approxi- 
mate the  propagation  of  the  scattered  waves  from  the  en- 
sonified  region  in  terms  of  the  free  field  Greens  function: 

ikR2    .£   £i 
G(x2xi)s  ^—  e~1Xs   x  (47) 

R2 

where 

*2 

We  further  assume  that  each  surface  element  da1  acts  as  a 
small  reflector,  and  that  the  response  of  da1  to  the  inci- 
dent wave  is  that  of  a  "local  reaction,"  i.e.,  independent 
of  any  other  part  of  the  ensonified  area  A. 

Then,  the  Kirchhoff's  approximation  where  it  is  assumed 
that  p  and  3p/3n  vanish  everywhere  on  the  surface  except  at 
the  ensonified  area  and  that  the  values  of  p  and  9p/9n  are 
proportional  to  the  incident  wave,  allow  the  scattered 
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"target"  strength  for  the  ensonified  area  to  be  approxi- 
mated as  : 

P  (x1)  !  =  R  PMU)  (x1)  (48) 

A       w 

1P_  !  =  R^P(i) 


an1  A      an1 
where 

p  is  the  locally  reflected  wave. 
R  =  reflection  coefficient. 


p^^os  9-pc  cos  91  m       c    _   c l 


p^cos  9+pc  cos  91   s       sine1 
Finally,  assuming  a  Gaussian  illumination  function: 

D  =  e  X2    Y2  (49) 

where 

X  and  Y  are  the  effective  dimensions  of  the  illuminated 

area. 
The  scattering  integral  based  on  the  Kirchhoff's  approxi- 
mation can  be  expressed  as: 

ik(Ri+R2)  f[      2i(ax1+6y1) 

Pw    (x^  =  -  2ttRxR2 Rf  (8x8,9,)  ^  De 

e2iC(x  y  )  dxidyi  (50) 


where 


a  =  j  (sin  9i  -  sin  92  cos  03) 
8  =-  5-  sin  92  sin  93 
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Y  =  -^-(cos  9!  +  cos  82) 

^(x^1)  =  surface  displacement  function. 


The  above  scattering  integral  is  then  applied  to  a  randomly 
rough  surface  where  the  surface  displacement  function  £  is 
a  random  variable  assumed  to  be  represented  by  a  Gaussian 
PDF  expressed  as 

W(;)  =  -±—   e"c2/2a2  (51) 

a/2T 

with  zero  mean  and  variance  a2. 

Also,  assume  that  the  surface  is  slowly  varying  so  that 

the  signal  reflects  from  an  essentially  stationary  surface 

and  that  the  succession  of  received  scattered  signals  p 

n 

are  assumed  to  form  a  satistically  independent  set  of 
sample  functions  from  which  sequence  N  the  first  and  second 
moment  of  the  field  can  be  investigated.   Doing  so  Clay  and 
Tolstoy  [Ref.  15]  found  that  the  mean  reflected  signal  can 
be  expressed  as 


<P  (s)(x2)>C  =  P„e"2Y2a2  (52) 


where 


Po  ■  PJS>  (x*>  I  ==0 

is  the  signal  reflected  by  a  mirror-like  surface, 

other  factors  being  the  same. 
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It  is  seen  from  the  above  that 

1.  For  a+0,  all  displacements  £  have  zero  proba- 
bility and  the  mean  signal  tends  to  p  .   Further- 
more, all  elements  contribute  to  the  scattering 
coherently. 

2.  For  a>>k,  all  displacements  5  are  equally  probable. 
There  are  large  phase  shifts  between  contributions 
from  different  surface  elements  and  they  tend  to  cancel 
each  other  and  the  scattering  radiation  is  incoherent. 


T 
<s2(t)  -  s  /"  s2(t)dt  =  <pp*>   -<pp*>         (53) 


The  second  moment  is  defined  as: 

1   T 

I  /   s2(t)dt  =  <pp*, 

0  L 

where 

1  N 
s2(t)  =  ±[l   Pn2(t)  -  N<p(t)>N2]  (54) 

In  this  expression  the  operation  of  squaring  the  signal  has 
to  be  considered.   However,  each  surface  element  has  a 
different  ;  for  a  random  surface.   Thus,  the  probability 
of  finding  element  dx *  dy'  with   £  1 ,  and  element  £ 2  with 
dx"  dy"  is  expressed  in  the  bivariate  distribution  function 
assumed  to  be  Gaussian  and  of  the  form: 


WttiCa)- re[ (Ci2  +  Ci2-2CiC2^)]  (55) 

27ra2(l-'j>2P    2(l-^2)a2 


where 


•M£,n)  =  —  <Ci(x1y1t)  C2(x"y"t)>  (56) 

a2 
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is  the  cross-correlation  function  characterizing 

the  surface  shape. 
Changing  to  polar  coordinate  leads  to  Eq.  6.51  of  Ref.  15 
where : 

<s2>a  /  D  Jo(2*r)[e-4Y2a2(1-*>-e-4^°2]  rdr       (57) 

Here  JQ  is  the  Bessel  function  of  zero  order  and  x  is  the 
transformation  parameter  given  as 

<=    a 


cos  9 
As  seen  from  Eq.  (57) 

1.  For  rough  surface,  yzoz  is  large  and  the  second 
term  in  the  bracket,  the  coherent  part,  is 
negligible. 

2.  For  smooth  surface,  y2o2is  zero  and  the  whole 
bracket  is  zero. 

Since  Eq.  (57)  cannot  be  integrated  directly,  Clay  and 
Tolstoy  [Ref.  15]  consider  it  for  small  and  large  ya 
separately. 

Thus,  concentrating  on  the  high  frequency  limit,  Clay 
and  Tolstoy  assumed  y2o2>>l   and  showed  that: 

-The  coherent  component  is  negligible. 

-<s2>  =<pp*>,  as  the  means  tends  to  zero, 

00 

<s2>=<pp*>ajjo(2Jcr)e"4Y2a2(1",|j)rdr;  y2a2>>l    (58) 
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Equation  (58)  consists  of  the  product  of  an  oscillatory 
function  and  an  exponential  function.   Because  of  the  Bessel 
function,  the  main  contribution  to  the  integral  is  near  r=0. 
Near  r=0  the  phase  changes  slowly  and  the  expression  can  be 
evaluated  by  the  method  of  stationary  phase.   Thus,  the 
expression  of  ty  about  r=0  is  given  as: 

ty    z    1+iF  (0)^-  (59) 

Furthermore,  Clay  and  Tolstoy  show  that  ^"(0)  can  be  related 
to  the  characteristics  of  the  surface  as: 


^■.(O)  |   =  A-  <C12>  (60) 

a 

Finally,  Clay  and  Tolstoy  show  that  the  scattering  signal 
can  be  expressed  as: 

4yza2>>l 
<s2>  =  <pp*>  =  <PlPl*>   -£-  Shf;  (61) 

R22  {2y2<^l2>)>>l_ 

A  =  ensonified  area 

<PiP?>:   the  expected  average  value  of  px2,  where 

Pi  is  the  incoming  pressure  to  the 

luminated  area. 


where 


S.  -  -  scattering  function 
hr 

x 


fZ(9)* _e  (2Y2<C12>) 


2tt(cos  9i  +  cos  92)  2<Ci2> 
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l+cos0 icos62-sin0 isin92cos0  3 
f  (0)  =  f  (0!0203)  = 1 

COS0  1+COS0  2 


R  _  P  C  COS0  -  QC  COS01      c         c1 

p  c  cos0 +  pc  cos©1   sin0    sin©1 


Y  =  -j    (cos  0!  +  cos  02) 


K  =  _iL 


COS 


a  -  tz    (sin  0i-  sin  02  cos  03) 


k  =  21 

K  X 

<Ci2>  =  (3  +  5.12W)  x  10"3;  W  =  wind  speed  in  m/s. 
Equation  (61)  is  valid  for: 
4y2a2>>l 
(2y2<Cl2>) >> 


In  summary,  the  reflection  of  high  frequency  signals  yields 
scattered  radiation  which  is  incoherent.   Furthermore,  as 
pointed  out  by  Clay  and  Tolstoy  in  [Ref .  15]  although  the 
radiation  is  primarily  scattered  in  the  specular  direction, 
parts  are  scattered  in  all  directions.   As  seen  from 
Eq.  (61)  the  scattering  function  S,  f  is  primarily  dependent 
on  the  mean  square  slope  of  the  surface  <£12>  and  neither 
the  mean  square  wave  height  a2  nor  the  correlation  distance 
are  important.   Furthermore,  it  is  noticeable  that  in  the 
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high  frequency  limit  the  scattering  function  is  independent 
of  the  frequency  since 


e  -  +1   when  r2-*» 

(2Y2<C12>) 


C.   GEOMETRICAL  SHADOWING 

The  phenomena  of  shadowing  of  some  surface  areas  by 
others  has  to  be  considered  either  when  the  surface  irreg- 
ularities are  large  with  respect  to  the  incident  wavelength 
or  when  the  grazing  angle  is  small. 

The  few  papers  devoted  to  this  subject  are  aimed  mainly 
towards  calculation  of  a  "Shadowing  function"  based  on 
the  statistics  of  the  surface. 

An  explicit  method ,  geometrical  shadowing,  has  been 
introduced  by  P.  Bechmann  [Ref.  17]  where  the  shadowing 
function  SO)  is  the  probability  that  the  point  ?  (Fig.  25) 
is  illuminated. 

00 

S(9)  =  exp  [-/q(x)dx]  <62) 

o 

where : 

q(x)  is  the  probability  that  C  is  shaded  by  z,   in 
the  interval  (x,x+dx)  given  that  it  is  not  shaded 
in  (0,x) . 
This  calculation  of  SO)  only  considers  the  elevation  of 
the  surface  observation  point.   However,  the  slope  also 
plays  a  role  in  that  if  its  value  exceeds  cote  the  point 
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will  certainly  be  shaded.   Thus,  R.  Wagner,  in  Ref.  18, 
incorporated  both  5  and  c1  using  P.  Bechmann's  method  and 
found  that  the  conditional  probability  that  a  point  on  the 
surface  is  illuminated,  given  that  it  has  height  5  and 

slope  5  ,  can  be  expressed  as: 

00 

SOiUi/Ci1)    =   exp[-  J  q(x)dx]       u(cot0    -    ? 1 )  (63) 

0 

where 

u:   is  the  unit  step  function. 

q(x):   is  the  conditional  probability  that  5  is 

shadowed  in  the  interval  (x,x+dx)  given  that 
it  is  not  shadowed  in  (0,x). 
The  function  q(x)  cannot  be  calculated  exactly.   Thus, 
R.  Wagner  made  the  approximation  that,  for  all  x,  the 
probability  that  5  crosses  the  ray  in  dx  is  independent  of 
the  values  of  C  and  c1  at  x=0.   In  the  above,  no  mention 
has  been  made  of  the  direction  of  observation.   However, 
in  this  respect,  R.  Wagner  [Ref.  18]  points  out  that  in 
the  high  frequency  limit  only  those  portions  of  the 
surface  which  are  illuminated  simultaneously  by  rays  in 
the  direction  of  incidence  and  observation  can  contribute 
to  the  observed  scattered  power.   For  this  condition, 
R.  Wagner  introduced  both  directions  as  independent 
variables  in  the  so-called  bistatic  shadowing  probabilities 
where  he  defines: 
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1.  S(9i92|Si,£i')  to  be  the  conditional  probability 
that  the  surface  will  not  cross  the  incoming  ray 
(Ray  1)  or  the  outgoing  ray  (Ray  2)  anywhere, 
given  that  both  rays  pass  through  an  arbitrary 

point  or  the  surface  with  displacement  £  and  slope  z, l  . 

2.  S  (9i  |  9 2 / C  i  / s  i  ' )  to  be  the  conditional  probability 
that  the  surface  does  not  cross  Ray  1,  given  that 
it  does  not  cross  Ray  2  and  that  both  rays  pass 
through  the  point  z,   having  slope  c '  . 

Thus,  the  conditional  shadowing  function  can  be  expressed 


as 


s(e1,e2|s1,c1')  -  scejez/CifCi')  s(e2|'ei,Ci')  (64) 

The  shadowing  function  is  then  obtained  by  averaging 
over  all  possible  heights  and  slopes 


00 

(91,9a)  =  //s(91/92U^')  W(5,£')d?  &V 


(65) 


Here  w(c,£')  is  the  bivariate  PDF  of  the  surface  height 
and  slope,  assumed  to  be  Gaussian 


*<?,?')  =  £  <*0I*0"I>"*  exp{-i-  -rffrTT'      ("> 


1 

0 

where  '>o  =a2  and  0«  ■  are  the  values  at  t=0  of  the  correla- 
tion function  and  its  second  derivatives,  respectively. 


For  the  region  0<9<J  where  the  probability  of  crossing 
one  ray  is  assumed  independent  of  that  of  crossing  the 
other,  R.  Wagner  found  that  the  bistatic  shadowing  function 
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could  be  expressed  as: 


{l-exp[-2(B!+B2) ]}x{erf  Vj+erf  v2 } 

3(61.62)  = — — t  an  \ 

4(Bi+B2)  (67) 


where 


exp  (-v.  2)-/rF  v,  erfc  v, 
4tt  v. 


1       1       1 
Bi ; ;   i  =  1,2 


1 


V-  =  =  —±-    ;  1  =  1,2 

2a2^  "|    2Z2 


as  we  from  Eq.  (60)  have  that  o2\\p"{0) 
=    <c,lz>   =  Z2 


r}^   =   cot  6;  i  =  1,2 

and  noting  that 

2 

Error  function  erf(x)  =  — 

/¥ 

/   -t2 
f    e  ^  dt 

0 

Complementary  error  function 

2 

erfc (x)  =  — 

7T 

J 

X 

The  shadowing  function  S(9i,62)  is,  in  short,  the  fraction 
of  the  surface  still  illuminated.   As  seen  from  Eq.  (60) , 
the  scattered  field,  in  the  high  frequency  case,  is  pro- 
portional to  the  illuminated  area.   Hence,  the  shadowing 
effect  of  a  rough  surface  can  be  introduced  by  multiplying 
the  ensonified  area  A  by  the  shadowing  function  S(01,e2). 


60 


D.   ESTIMATING  THE  SURFACE  SCATTERING  EFFECT 

The  following  estimations  are  based  on  calculations 
in  the  specular  direction,  which  is,  as  pointed  out 
earlier,  expected  to  give  the  maximum  supplementary 
scattering  effect.   Hence,  in  the  specular  direction  where 
81  =  3  -  =  5  and  9 3  =  0  ,  the  scattering  function  S,  ~ 

reduces  to  the  following  expression: 

2 

f2(9)  R2  K- 


as 


S.  -  =  *-    ^'    - e   (2y2<ci2>)       (68) 

hf    2  (cos0i+cos92) 2<c12> 

s   =   R2 

hf    8tt<C12> 

f(6)  =  cos9 

Y  =  k.cos9 

a  =  0 

x  =  -£L-=  0 
cos9 

As  pointed  out  in  the  previous  paragraph,  the  reflection 

of  very  high  frequency  signals  by  the  sea  surface  yields 

scattered  radiation  that  is  incoherent  under  the  assumption 

that 

4yza2>>l 
(2Y2<C12>)»iT 
Before  we  launch  into  the  calculations,  we  will  verify 
these  criteria  for  the  frequency  range  of  interest:   60  and 

30  kHz. 
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Utilizing 

Y2  =  k2  cos29  a  =  8.1  x  10"3 

?     aW1* 

a   =  ~  6  =  0.74 

4Sg2 

<C12>  =  (3  +  5.12W)  x  10"3    g  =  9.81  m/s2 
and  assuming  a  windspeed  of  10  m/s  (SS3)  yields 
a2  =  0.30  m2  ;  a  =  0.55  m 

<£12>  =  5.42  x  10~2 
For  60  kHz (A  =  .025  m)  and  low  grazing  angles,  e.g.,  9  =  85( 
4y2a2  =  5.75  x  102  »1 

(2Y2<C12>)  =  5.2  x  101  >>  i 

R2 

and  the  first  order  5.2>>—  ,  R>0.5m. 

R2 

Thus,  for  the  60  kHz  case,  the  criteria  are  fulfilled. 

For  30  kHz  (X=.05m)  and  9=85°. 

2„  2    _    1     A    A  ..    T  n2. 


4cTy   =  1.44  x  10"»1 
(2y2<C12>)  =  13>>  i 


R 
again  the  first  order   R>1.0  m. 

Thus,  also  for  the  30  kHz  case  the  criteria  are  fulfilled. 
Similarly,  for  the  shadowing  function  in  the  specular 
direction 

0<9i  =  9  2  =  9  <  tt/2    and  9  3  =  0 
we  obtain  the  following  simplified  expression 
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S(0)  =  S(81,82)  =  [l-exp(-4B)]erfv 

4B 


(69) 


as 

Vl  =  v2  =  v  = 


(2a2|i|;o"|^      (2>2)h 


n  i  =  n 2  =  n  =  cot9 

B   =  B   =  B  =  [exP(~v'2)~>/^  v  erfc  v] 

4/tt  v 

In  summary,  for  specular  scattering  at  the  high  frequency 

limit,  the  expected  average  value  of  p2,  where  p  is  the 

pressure  field  at  the  receiver,  is  then  obtained  from  the 

following  simplified  expressions: 

<s2>  =  <pp*>  hSIZL   shf  (70) 

R2 


2»     _ 


<S'>    =    <pp* 


fc      AS(9)R: 


2 


1  1 


8ttR2    <X,      > 


where 


_...    _    l-exp(-4B)erfv 
S(9) 4B 


v  =   ^^^ 


(2ZM 


=   cot9 


2  \  *2 


[exp  (-v2  )  -/tTv  erf  c  v] 


B   -  -5 

4(tt)  2   v 
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<-12>  =  i2    =  (3  +  5.12W)  x  10" J  ;  W  =  windspeed  in 
m/s. 

R   p^cose  -  p  c  cose1  m     _c_     c 1 


P^cose  +  p  c  cosG1    siny    sine1 

From  Ref.  19  we  use  the  following  air/sea  water  inter- 
face data: 

1.  For  air  c'  =  343  m/s  and  p^1  =  415  Rayls. 

2.  For  sea  water  c  =  1500  m/s  and  pc  =  1.54  x  10   Rayls. 
The  estimation  of  the  illuminated  area  A  for  specular 

scattering  where  Q i  =  9  2  =  0  and  e3  =  0,  is  based  on  the 
geometry  illustrated  in  Fig.  24b.   The  illuminated  area 
is  given  by 

A  =  uab  (71) 

Assuming  the  following  data  to  be  known 

h  =  depth  of  the  source 

h  =  depth  of  the  receiver 

r 

A4>  =  the  half  beam  width  of  the  directional  receiver, 
Both  a  and  b  of  Eq.  (71)  can  be  calculated  in  terms  of  the 
detection  range  R_  as  follows: 


x2  = 


h  +h 
s   r 


(72) 


h  _x    RD 

xi  =  Rp  -  x2  =  RD(l-H-fh-)  ;  6  =  tan   (h-Th-) 

s   r  s   r 
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Ri 

X 

1 

sine 

R2 

X2 

To  a  good  approximation  when  9  =  tt/2  ,  we  have 

hr 

a  =  Xi  =  RD(1  "  h-Th-)  <73> 

s  r 

b  =  R  sin("*>  ■  (h  +£  )sine  sin(A*> 

s   r 

The  expression  for  <s2>  is  then  introduced  as  a  supple- 
ment to  the  direct  path  to  the  receiver  in  the  following 
way: 

By  utilizing  the  relationship 

1  =  h  <74) 

pc 
for  the  intensity,  the  scattering  intensity  at  the  receiver 
and  the  intensity  at  the  ensonified  area  are,  respectively, 

I   =  IE  _  <sz> 
s    pc     c 

p2   <P1P1*> 

T     _   _   _ 

1     PC       PC 

Equation  (70)  can  thus  be  written  as 
t     t   AS  ( 9 )  R2 

S       *         2    l  •} 

s       8ttR2  <c12> 


_  ._      AS  (9)  R2 


S  87rR22<<;12> 
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lOlog    I   /I,    =    lOlog    /^ref   =   K 

Ii/Iref  c 


ILs    =    IL(Rl)    =    Kq  (75) 

We    then   have   to   determine    IL(Rj) 

TL    =    SPL(l)     -    SPLtRj)    =    201ogR1  (76) 

SPLCRJ    =    ILCRJ    =    lOlogl! 

SPL(l)     =    SL 
yielding 

IL(Rt)  =  SL  -  201ogRi  =  K:  (77) 

Then 

IL     .  -   K   +  K,  (78) 

suppl     o     l  v   ' 

is  the  supplement  to  the  direct  path,  and 

IL , •    ,  =  SL  -  201ogR^  =  K, 
direct  '  D     3 

x       /t        j.  ■  t      direct    „ 
rdirect/3Cref  =  antl  log  To =  K<* 

IL 

I  i/I      jt   =   anti    log   -rr-  =   K5 

suppl7     ref  ^10  5 

The   total    intensity  of   the   receiver   is   thus 

t        /I  =   direct  +   ^uppl   m 

Tot7    ref  I      -  4  5  6 

ref 


ILTQt  =  101ogK6  =  K7  (79) 
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Thus,  the  effect  of  the  randomly  rough  surface  compared 
to  the  idealized  free-field  condition  can  be  expressed  as 

AIL  =  IL_  .  -  IL, .    .  =  K8  (80) 

Tot     direct     8 

A  calculator  program  on  a  Texas  Instrument  5  9  (later 

called  TI  59)   was  developed  to  perform  these  calculations. 

A  block  diagram  of  the  program  is  outlined  in  Fig.  (26) 

and  the  programs  steps  together  with  a  detailed  description 

is  given  in  Appendix  B.   The  calculations  are  based  on  the 

following  fixed  data 

h  =  2  m 

s 

h   =  6  m 

r 

A<J>  =  10° 

W   =  10  m/s  (SS3) 
Then,  varying  the  detection  range  from  R  =  2000  m  to  100  m 
gives  the  difference  between  ILT   and  ILdj_rect  plotted  in 
Fig.  (27). 

E.   CONCLUSION  AND  DISCUSSION 

As  seen  from  Fig.  (27)  the  supplementary  effect  of 
the  scattering  from  a  rough  surface  in  the  high  frequency 
case  is  negligible  compared  to  the  direct  path. 

In  saying  so,  it  also  should  be  pointed  out  that  the 
Helmholtz-Kirchhoff 's  approach  may  be  limited  as  it  does 
not  take  into  account  the  diffraction  effects  from  crests 
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and  throughs  of  the  ensonified  area,  an  effect  which 
becomes  increasingly  important  at  low  grazing  angles, 
high  frequency  and  when  the  rough  surface  is  a  super- 
position of  swell  and  capillary  waves. 
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VII.   THE  EFFECT  OF  SCATTERING  AND  ABSORPTION  FROM  THE  SUB- 
SURFACE OCEAN  LAYER 

A.   GENERAL  SCATTERING  THEORY 

Generally  when  a  region  (volume)  scatters  sound,  some 
of  the  energy  carried  by  the  incident  wave  is  dispersed. 

The  energy  lost  by  the  incident  wave  may  be  absorbed 
by  the  scatterers  or  it  may  be  simply  deflected  from  its 
original  course.   The  amount  of  energy  lost  per  second  by 
the  incident  wave  divided  by  the  incident  wave's  intensity 
is  called  the  total  cross  section  a   of  the  region  and  is 
the  sum  of  the  absorption  and  scattering  cross  sections 

ns+Ea 


where 


II   =  scattered  power 


II   =  absorbed  power 
a 


The  existence  of  gas  bubbles  in  the  subsurface  ocean 
layer  modifies  the  forward  scattering  in  the  following 
two  major  ways: 

1.   The  bubbles  can  resonate.   When  the  bubbles  are 

excited  at  a  frequency  near  its  natural  frequency, 
it  very  efficiently  absorbs  and  scatters  the  inci- 
dent wave.   At  resonance,  the  scattering  and 
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absorption  cross  section  of  a  typical  bubble  at  sea 
is  of  the  order  10   times  its  geometrical  cross 
section. 

2.   The  bubbles  change  the  effective  compressibility  of 
the  water  and  cause  the  speed  of  sound  to  be  a 
function  of  frequency,  i.e.,  the  medium  is  dispersive. 

We  will  investigate  and  discuss  these  effects  by  separately 
estimating: 

1.  The  attenuation  due  to  the  bubbles  and 

2.  the  refraction  by  bubbles. 

B.   ABSORPTION  MODEL 

The  choice  of  model  for  the  subsurface  ocean  layer  depends 
on  whether  the  medium  has  a  teneous  or  a  dense  distribution 
of  scatterers.   When  the  bubble  density  is  teneous,  both 
"single  scattering"  and  "first  order  multiple  scattering" 
approximation  solutions  are  applicable. 

On  the  other  extreme,  when  the  bubble  density  is  high, 
the  so-called  "diffusion"  approximation  can  be  used. 
Between  these  two  extremes,  multiple  scattering  effects  are 
important . 

The  multiple  scattering  theory,  which  in  the  limit  also 
contains  the  first  order  approximation,  will  be  used  to 
estimate  effects  of  attenuation  due  to  bubbles,  on  the 
propagation  of  propeller  noise  from  the  target  to  the 
torpedo. 
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The  geometry  of  the  propagation  model  is  illustrated 
in  Fig.  (28)  where  it  is  assumed  that  a  plane  wave  is 
incident  on  a  semi-infinite  (disregarding  the  sea  surface) 
slab  of  thickness  x  containing  a  number  of  randomly  dis- 
tributed bubbles.   The  plane  wave  approximation  is  valid 
if  the  incident  sound  has  a  wavelength  X   much  greater 
than  the  bubble-radius  a 

ka<<l 
where 

k=aj/c=2T/X 
The  receiver  is  located  outside  the  slab  and  the  beam 
pattern  of  the  receiver  is  represented  by  the  solid  angle 

Q    . 

r 

We  are  interested  in  the  estimating  of  the  total  power 
received,  taking  into  account  the  multiple  scattering 
process  in  the  inhomogeneous  slab  as  well  as  the  beam 
pattern  of  the  receiver. 

The  mathematical  formulation  of  this  problem  is  based 
on  Twersky's  theory  of  multiple  scattering.   Since  the 
theory  is  presented  in  Ref.  20  only,  the  basic  formulation, 
major  assumptions,  and  the  end  results  will  be  presented 
here . 

The  total  intensity  is  the  average  of  the  square  of 
the  magnitude  of  the  total  field: 
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<|*a|2>    =    <|<^a>    +  ^fa|2>   =    |<^>|2    +    <|*   *|a>       (81) 


where 


**  =  #  *  +  £   u  a 
1     8-1   S 

the  scalar  field  at  the  receiver  location  r  , 

a' 

see  Fig.  (29),  is  the  sum  of  the  incident  wave  d> . 
and  the  contribution  from  all  N  scatterers. 
|<iji  >|2  is  the  coherent  intensity  based  on  the 
average  field  <ij;a>.   <|ipfa|2>   is  the  incoherent 
intensity  based  on  the  fluctuating  field  ij;  a 
In  Twersky's  theory,  the  multiple  scattering  process  is 

described  by  the  following  set  of  integral  equations  which 

Eq.  (81)  must  satisfy: 

<|*a|2>  =  |<*a>|2  +   |v  a|2  <|*s|*>p(r"  )drQ    (82) 


where 


vsa=  Us3+  "t3  vst  P«?t>«* 


is  an  operator  representing  all  the  scattering  processes 
from  s  to  a.   (See  Fig.  (28).)   It  should  be  noted  that 
Twersky's  theory  includes  all  the  multiply  scattered  waves 
that  involve  chains  of  successive  scattering  going  through 
different  scatterers.   (See  Fig.  (29a).)   However,  the  theory 
neglects  the  terms  which  include  an  individual 
scatterer  more  than  once,  as  illustrated  by  Fig.  (29b) . 
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Thus  Twersky's  theory  is  expected  to  give  good  results 
when  back  scattering  is  insignificant  compared  to  the 
scattering  in  other  directions. 

As  typical  for  most  integral  equations,  Twersky  being 
no  exception,  detailed  complete  solutions  are  difficult 
to  obtain.  However,  Twersky  gave  an  approximate  solution 
to  Eq.  (81)  and  Eq.  (82) ,  which  according  to  Ref .  20  have 
been  found  to  agree  reasonably  well  with  experimental  data 
This  solution  is  based  on  the  following  main  assumptions: 
-Backscattering  is  assumed  to  be  small  compared  to 

scattering  in  other  directions. 
-Scattering  is  mostly  concentrated  in  the  forward 
direction.   This  is  reasonable  based  on  the  assumed 
random  distribution  of  the  bubbles;  i.e.,  no  rein- 
forcement of  the  radiation  pattern  occurs  except  in 
the  direction  of  the  incident  wave. 

-The  angle  9    is  small,  i.e.,  9    =  0  . 
3     as  as 

This  leads  to  the  following  expression  for  the  total 
intensity  at  the  receiver: 

<\\ba\2>   m   T  -   exp(-pa  x)  exp(-pa  x)  (83 

+q[l-exp(-pasx) ] 


where 


a   =  absorption  cross  section 
a 

a      =  scattering  cross  section 

s 
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Q  '      s 
q=  r 


,.r|f  I  2dft 

4tt  '  '    s 

where 

f  =  the  amplitude  function 

q  =  the  fraction  of  total  scattered  power  collected  by 

the  receiver  as  illustrated  in  Fig.  30  and 
p  =  is  the  bubble  density,  i.e.,  the  number  of 

scatterers  per  unit  volume. 
For  small  values  of  px  we  see  from  Eq.  (82)  that  the 
coherent  part  dominates: 

InT  =  ~{o+o) px  (84) 

In  this  limit  the  multiple  scattering  result  is  equiva- 
lent to  that  obtained  from  single  scattering  considerations 

For  large  values  of  px  which  corresponds  to  very  dense 

or  very  wide  slab  of  scatterers,  the  incoherent  intensity 

dominates 

InT  =  lnq  -  a  px  (85) 

a 

In  this  latter  case,  it  is  notable  that  when  9.   =2tt, 


the  receiver  collects  almost  all  the  scattered  power 

InT  ;  -  a  px 
a 


(86) 


The  first  case,  representing  the  situation  for  teneous 
density  of  scatterers  and/or  narrow  beam  pattern  of  the 
receiver,  gives  a  good  approximation  to  the  situation  of 
interest  in  the  thesis.   It  also  represents  the  case  for 
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which  no  scattered  power  is  received.   This  will  be 
approximately  true  for  a  narrow-beam  width  receiver.   We 
will,  therefore,  investigate  the  coherent  intensity  first. 

Also,  the  incoherent  case,  as  represented  by  Eq.  (86) 
will  be  investigated,  where  only  losses  due  to  absorption 
are  incorporated. 

C.   THE  COHERENT  INTENSITY  CASE 

For  the  coherent  case  the  intensity  level  after  the 
incident  wave  has  traversed  a  distance  x  is 

Ix  =  I   exp[-(aa+as)px]  (87) 


where 


I   =  incident  plane  wave  intensity. 
P 


The  change  in  intensity  over  the  distance  x  is 

Ix   I  exp  [-(aa+ac) px] 


AIL  =  lOlog^  =  "P  —  ^  x  a   S" 


=  lOlogexp  [-(a  +a  )px] 

a.         s 

The  excess  attenuation  per  unit  distance  due  to  bubbles 
is  thus 


a  =  - 


AIL      r  x  a   s 


-  (a   +a    )  px  log  e 


=  -  [ — - ] 


=  (a  +a  )  p log  e 
a  s 

a  =  4.34  a   p,  in  dB/m.  (89) 

e 

However,  this  only  takes  into  account  bubbles  of  one 
size.   In  a  bubbly  medium  there  is  a  spectrum  of  radii. 
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The  probability  density  function  for  finding  a  bubble  size 
between  radii  a  and  a+da  is 
W(a)  =  2i*L 


P 


(90) 


where 


/  W(a)da  =  1 

o 

n(a)da  is  the  number  of  bubbles  per  unit  volume 
having  radii  between  a  and  a+da.   It  is  common  to 
use  da=l  ym. 


n(a)da  is  the  total  number  of  bubbles  per 
o 
unit  volume. 

As  the  extinction  cross  section  also  is  a  function  of 

the  radius    (See  Eq.  (93).),  the  absorption   due  to 

bubbles  is  obtained  by  integrating  Eq.  (89)  over  all 

possible  radii 

oo 

a  =  4.34  J    a  (a)n(a)da  (91) 

o 

To  calculate  the  absorption  coefficient,  the  extinction 
cross  section  a   must  be  derived  from  the  general  bubble 
dynamic  relationship.   This  is  done  in  detail  by  C.  Clay 
and  H.  Medwin  [Ref.  21]  from  which  the  following  results 
are  taken. 
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The  scattering  cross  section 

n 

°s   =   T-  (92) 

p 

4ira2 


where 


[(fr/f)2-i]2+s 


f   =  resonance  frequency  =  ^i—  ( — ) 

f  =  f   =  operating  frequency 


2yb3P 


A,  h 


6  =  damping  constant  =6   +  6,  +  5 

r    t    v 

=  ka  +  <<*)  <^)2  +  -St- 

b   f      PAwa2 

a   =  bubble  radii 


—    3h-l)  r2i ( sin-  h  X  +  sin  X)  -2  (cos  h  X  -  cos  X) -, 

X2 (cos  h  X  -  cos  X)  +  3(7-1) (sin  h  X  -  sin  X 


b  -2 


3w  C    , 
X   =  a(   gg  P<3)h 

g 

k   =  thermal  conductivity  of  gas 

2t 

p   =  density  of  gas  =  pgA[l+^-§y]  (1+O.lz) 

p   =  density  of  gas  at  sea  level 

t  =  surface  tension 

P,  =  1.013  x  106(l+0.1z) 
A 

z  =  bubble  depth  in  m. 

C   =  specific  heat  of  constant  pressure  of  gas 

pg 
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U  =  shear  viscosity  of  water 
Y  =  7/5,  for  diatomic  gas 
DA  =  density  for  sea  water 


and 


b  =  [i+(|)^]-i  [x  +  3y-l   sin  h  X  -  sin  X] 

X   cos  h  X  -  cos  X 


8  "  *  +  ip  «  -  37b» 

Furthermore : 

0=0+0=        47ra2(6/ka) (93) 

S    a   S  [(fr/f)2  -l]2+62 


a      =o-a  (94) 

a    e   s  v   ' 


A  detailed  computer  program,  as  outlined  in  Appendix  B, 

was  developed  for  the  TI  59  to  handle  the  derivation  of 

a   ,    0    ,    and  a      based  on  an  assumed  receiver  depth  of  z  =  6  m. 
s   e       a 

For  the  60  kHz  case,  both  the  extinction  cross  section 

j   and  the  absorption  cross  section  a   are  given  in  Fig.  31 
e  a 

as  a  function  of  bubble  radius  a. 

Similarly,  Fig.  32  gives  a   and  a   for  the  30  kHz  case. 
■* '    3     3      e      a 

Superimposed  on  these  figures  are  the  curves  for  n(a)da  as 
calculated  from  the  following:   Figures  5  -  7  of  Ref.  2 
give  the  resonant  bubble  densities  in  a  1  urn  band  as  a 
function  of  depth  and  with  the  wind  speed  as  parameter  for 
the  three  discrete  frequencies  12  kHz,  38  kHz  and  120  kHz. 
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Based  on  these  data,  Fig.  33  shows  the  interpolated  bubble 
density  as  a  function  of  resonance  frequency  for  sea  state 
2 ,  3 ,  and  6 . 

Furthermore,  both  A.  Lovik  [Ref.  2]  and  H.  Medwin  [Ref.  3] 
found  that  the  bubble  density  function  n(a)  decreases  with 

increasing  bubble  radii  as 

-x 

naa 

where  H.  Medwin  [Ref.  3]  found  the  power  law  of: 

x  =  4  for  a<50-80  pm. 

x  =  2  for  a>50-80  ym. 
and  A.  Lovik  [Ref.  2]  found  the  power  law  of: 

x  =  4.2  between  38  kHz  and  120  kHz 

x  =  2.6  between  12  kHz  and  38  kHz 
which,  averaged  over  the  depth  interval,  corresponds  to  the 
bubble  radii  of  380  ym(12  kHz),  120  ym  (38  kHz)  and  49  ym 
(120  kHz) . 

As  suggested  by  A.  Lovik  [Ref.  2] ,  the  discrepancy 
between  the  two  observations  is  not  great  and  may  be  due 
to  the  few  measuring  frequencies  used  in  the  work  of 
A.  Lovik  . 

In  summary,  Fig.  33,  from  which  we  obtain  the  appro- 
priate resonant  bubble  density  in  a  1  pm  band  n(a_) 
together  with  the  power  law  naa   ,  comprise  the  full 
knowledge  of  n(a)da. 
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Performing  a  multiplication  of  0q   and  n(a)da,  we 
obtain  Fig.  34  and  Fig.  35  for  the  60  kHz  and  the  30  kHz 
case,  respectively. 

Finally,  the  integral 

/  ^e(a)n(a)da 

o 

was  evaluated  using  a  numerical  integration  based  on 
Simpson's  discrete  approximation  programmed  for  the  TI  59 
and  documented  in  Appendix  C. 

Based  on  the  above,  the  following  absorption  coeffi- 
cient for  the  coherent  case  is  obtained  for  f  =  60  kHz 

oo 

a  =  4.34  /  a  (a)n(a)da  =  4.34  (2.016  x  10"1) 

o    e 

=  8.75  x  10"1  dB/m. 


with 


a  *  0.8  8  dB/m  (95) 


n(aR)  =  1000 


n(a)  a  a-4 

z  =  6  m. 
Thus,  it  is  seen  that  the  attenuation  due  to  bubbles 
is  considerably  greater  than  the  "normal  attenuation" 
due  to  chemical  and  viscous  relaxation  processes  in  sea 
water,  which  for  the  60  kHz  case  is 

a  =  0.02  dB/m. 
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Thus,  the  total  absorption  coefficient  for  the  60  kHz 
case  is 

a  *  0.90  dB/m.  (96) 

For  f  =  30  kHz,  the  absorption  coefficient  due  to  bubbles 
in  the  coherent  case  is 

oo 

a  =  4.34  /  ae(a)n(a)da  =  4.34  (3.8375  x  lO-2) 
o 

=  1.665  x  10"1  dB/m. 
a  *  0.17  dB/m  (97) 

with 

n(aR)  =  20 
n(a)  a   a"2*6 
z  =  6  m. 
For  this  case,  the  power  law  dependence  of  n(a)  a  a 
gives 

a   =  1.655  x  10"1 
Thus,  the  difference  in  power  law  dependence  makes  no 
significant  difference  in  the  absorption  coefficient. 

The  absorption  coefficient  due  to  chemical  and  viscous 
relaxation  processes  is  at  30  kHz 
a  =  0.012  dB/m. 
The  total  absorption  coefficient  in  the  coherent 
case  is 

a  s  0.18  dB/m.  (98) 
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D.   THE  INCOHERENT  INTENSITY  CASE 

For  the  incoherent  case,  where  only  losses  due  to  ab- 
sorption are  included,  the  product  a   and  n(a)da  for  the 

a 

60  kHz  and  30  kHz  cases  are  given  in  Fig.  36  and  Fig.  37, 
respectively.   Performing  a  numerical  integration  based  on 
the  Simpson's  discrete  approximation  leads  to  the  follow- 
ing results: 

For  f  =  60  kHz,  the  absorption  coefficient  due  to 
bubbles  is 


=  4.34  /  a  (a)n(a)da  =  4.34  (1.66  x  10_1) 


,-1 

a 

o 
a  =  0.72  dB/m.  (99) 


with 


n(aR)  =  1000 

n (a)  a  a 
z  =  6  m. 
Adding  the  "normal  attenuation"  in  sea  water  for 
f  =  kHz,  yields  a  total  absorption  coefficient  of 

a  *  0.74  dB/m  (10°) 

For  f  =  30  kHz,  the  absorption  coefficient  due  to 

bubbles  is 

a  =  4.34  fa    (a)n(a)da  =  4.34  (2.975  x  10"2)    (101) 
^      a 

o 

a  z    0.13  dB/m 
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with 

n(aR)  =  20 

n  ( a )  a  a 

z  =  6  m. 

Adding  the  "normal  attenuation"  in  sea  water  for  f  =  30  kHz 
yields  a  total  absorption  coefficient  in  the  incoherent  case 
of 

a  ~  0.14  dB/m  (102) 

E.   SUMMARY  AND  DISCUSSION  OF  THE  BUBBLE  ATTENUATION 

As  seen  from  the  above,  the  attenuation  due  to  bubbles  in 

the  subsurface  ocean  layer  is  important  for  high  frequency 

and  high  Sea  State. 

Accounting  for  the  effect  of  bubbles  at  Sea  State  3  in 

summary  we  found  the  following  absorption  coefficients  in 

dB/m: 

f  =  60  kHz  f  =  30  kHz 

coherent   incoherent     coherent      incoherent 

a  =  0.88      a  =  0.72      a  =  0.17        a  =  0.13 

The  main  assumptions  were: 

-the  back  scattering  is  small  compared  to  that  in 

other  directions, 
-the  scattering  is  mostly  concentrated  in  the 
forward  direction. 
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-the  angle   9 ag  between  the  scatterer  and  the 

receiver  is  small,  i.e..  9    -  0°. 

as 

-the  depth  of  the  receiver  is  z  =  6  m. 
-the  sea  state  is  3. 
The  range  dependent  portion  of  the  passive  sonar 
Eq.  (1)  TL=-201ogR-aR  for  both  the  60  kHz  and  30  kHz 
scattering  results  can  now  be  compared  with  the  reference 
data  obtained  in  Section  V. 

For  f  =  60  kHz,  Fig.  22  gives  TL=-201ogR-aR  as  a 
function  of  R  with 

1.  a  =  0.021  dB/m,  the  "normal  attenuation"  due  to 
chemical  and  viscous  relaxation  processes. 

2.  a  =  0.9  dB/m,  the  total  attenuation  including  the 
effect  of  bubbles  in  the  coherent  case. 

3.  a  =  0.74  dB/m,  the  total  attenuation  including  the 
effect  of  bubbles  in  the  incoherent  case. 

For  the  same  example  as  in  the  reference  model,  Fig.  22 
yields  the  detection  ranges  for  60  kHz. 

R  =  1200  m  for  a  =  0.021  dB/m  (103) 

R  =  60  m  for  a  =  0.90  dB/m 

R  =  70  m  for  a  =  0.74  dB/m 
Not  surprisingly,  this  result  seems  to  exclude  the  possi- 
bility or  having  both  a  searching  and  attack  depth  near 
the  surface,  i.e.,  z  =  6  m,  for  a  torpedo  system  operating 
in  a  high  frequency  region,  f  =  60  kHz. 
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Similarly,  for  f  =  30  kHz,  Fig.  23  gives  TL=-201ogR-aR 
as  a  function  of  R  with 

1.  ct  =  0.012  dB/m,  the  "normal  attenuation"  in  sea 
water. 

2.  a   =  0.18  dB/m,  the  total  attenuation  including  the 
effect  of  bubbles  in  the  coherent  case. 

3.  a  =  0.14  dB/m,  the  total  attenuation  including  the 
effect  of  bubbles  in  the  incoherent  case. 

For  the  same  detection  example  as  in  the  reference  model, 
Fig.  23  yields  the  detection  ranges  for  30  kHz: 

R  =  2400  m  for  a  =  0.012  dB/m  (104) 

R  =  250  m  for   a  =  0.18  dB/m 

R  =  310  m  for   a  =  0.14  dB/m. 
Again,  the  bubbles  give  a  major  decrease  in  the  detection 
range.   A  detection  range  of  R  =  250  m  seems  marginally 
acceptable  as  the  turn  rate  requirement  for  a  pursuit 
homing  trajectory  may  become  excessive. 

The  above  results  are  summarized  in  Tables  III  and  IV 
for  the  60  kHz  and  30  kHz,  respectively. 

F.   THE  REFRACTION  BY  BUBBLES 

The  presence  of  bubbles  in  the  sea  water  affects  the 
speed  of  sound  (phase  speed)  primarily  because  of  the  change 
in  compressibility.   The  derivation  of  this  dispersive  rela- 
tionship on  the  sound  speed  has  been  done  by  H.  Medwin 
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[kef.  22].   H.  Medwin  showed  that  the  bubbles  with  resonant 
frequencies  greater  than  the  incident  frequency  decrease 
the  sound  speed,  while  bubbles  with  resonant  frequency  lower 
than  the  incident  wave  increase  the  sound  speed. 

Furthermore,  H.  Medwin  [Ref.  22]  predicts  the  sound  speed 
gradient  due  to  bubbles  as  a  function  of  depth  for  a  frequency 
range  and  wind  speed  compatible  with  our  domain  of  interest. 
He  found   the  gradients 

g  =  dc/dz   =  0.26  s"1  at  z  =  0  m. 

g  =  3c/9z  =  0.016  s"1  at  z  =  10  m. 

g  =  dc/dz   =  0.005  s"1  at  z  =  20  m. 
For  comparison,  the  sound  speed  gradient  due  to  pressure 
in  an  isothermal  water  is 

g  =  0.017  s"1 
This  shows  that  the  rays  in  the  top  10  m  are  strongly 
influenced  by  bubbles.   However,  with  respect  to  our  surface 
scattering  model  (Section  VI)  where  both  the  source  and  the 
receiver  are  situated  very  close  to  the  surface  and  where 
relatively  short  propagation  distances  are  encountered, 
this  refraction  phenomenon  is  assumed  to  have  negligible 
effect. 
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VIII.   THE  TURN  RATE  LIMITATION 

As  pointed  out  in  the  previous  section,  the  presence 
of  bubbles  near  the  surface  may  significantly  reduce  the 
range  at  which  the  target  can  be  detected. 

In  this  section ,  the  turn  rate  necessary  during 
pursuit  homing  at  the  previous  estimated  detection  ranges 
will  be  studied.   A  computational  procedure  will  be  used 
to  determine  the  range  of  angles  on  the  bow  (AOB)  of  the 
target  at  the  beginning  of  homing  which  lead  to  miss  on 
the  initial  attack,   A  trajectory  where  the  torpedo  velocity 
vector  always  is  directed  towards  the  instantaneous  target 
position  is  called  a  pursuit  homing  trajectory.   The 
derivation  of  the  pursuit  homing  trajectory  follows 
P.  van  Nostrand  [Ref.  23]  and  is  based  on  the  geometry  of 
Fig.  38,  where 

r  =  approach  angle,  i.e.,  angle  between  ship 
velocity  vector  and  the  line  of  sight. 

AOB  =  180  -  <J> 
The  equation  of  motion  is  obtained  by  taking  the 
component  along  r  and  the  normal  to  r,  yielding 

r  =  V      cos    <J>   -  VT  (105) 


ri   =   -V      sin   <$>  (106) 
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where 

r  =  range  rate 

• 

<J>  =  turn  rate 

Vs  =  target  speed 

VT  =  torpedo  speed. 

Dividing  Eq.  (105)  by  Eq.  (106)  yields: 

-    VT    1 

r  =  (V~  iTnT'~  cot  *>*   ;  r  ^° 
s 

and  defining 

V 


T 

s 


V"  =  P 


yields 

I =  (iil*  - cot  *>*  (107) 

Then,  introducing  the  initial  conditions:   r  ,  tj>  where 
r  :   initial  detection  range 

<p    :   initial  approaching  angle 

yields  _ 

(1+cos  t>) 

K  ■  ro— — -f^r  (109) 

(sin  <P    ) 

For  the  geometry  of  Fig.  38,  the  turning  rate  is  given  by 
Eq.  (105) . 
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V 
i   =   ~   sin  <p  (110) 

Substituting  Eq.  (108)  into  Eq.  (110)  yields: 

V  (1+cos  $)P 

i  =  -  — 5-5 —  (111) 

K(sin  $) 

From  Eq.  (108)  we  see  that  r-*0  as  $-»-0,  i.e.,   $  tends  to 

zero  as  the  torpedo  approaches  the  target  ship.   It  is 

further  of  interest  to  determine  the  limiting  value  of  the 

turning  rate  as  the  torpedo  approaches  the  target  for 

various  values  of  the  parameter  p.   This  is  done  by  taking 

the  derivative  of  j>  with  respect  to  <J>.   Thus,  from  Eq.  (109) 

we  get 

|$  =   -  -#(sin   4>)1-P    [l+cos<p]P    [2   cos    <j>  -  p]       (112) 
dtp  k 

For   l<p<2,    we    see   that   Eq.     (112)    is    zero   at 

cos    5>  =    P/2  (113) 

$>  =   cos"1  (P/2) 

with  an  associated  maxima 


V    (1+P/2)P[1-(P/2)2]1"P/2 
s 


(114) 


1  4>max '  K 

Furthermore,    the   turn   rate    for  p<2   at   impact   is    zero   as 
the    limit   of   Eq.     (109)    yields: 

lim   l=   _      s    lim(sin    (}))2~p[i+cos    cj>]P  -*0  (115) 

K 

$+0  <$>-+o 
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For  p=2  we  see  from  Dq.  (Ill)  that 

V                     4V 
lim  $=  -  ~  lim(l+  cos  4>)   = ^  (116) 

4i+0  $+0 

For  p>2  we  see  from  Eq.  (Ill)  that  there  is  no  zero 
value  of  d<J>/dcJ)   between  i  =  0   and  <p   =  180  ,  since  all  terms 
of  Eq.  (112)  are  nonzero  terms  for  0  <(J)<180o.   Furthermore, 
the  turn  rate  for  p>2  at  impact  is 


V  P 

It-     II      I      T  ,  •  (1+COS  <j>)   ■   _  /-!1-7\ 

!  lim  $|  =  |  -  -=-  lim ^-^1^°°  (117) 

4,-0  (sin  0) 

as  (sin  p}*00  since  (p-2)>0. 

Furthermore,  as  $-»-180  ,  we  get  from  Eq.  (112) 


•    V 

lim  ^  =  -#(2+p)  lim  (sin  $) 1"P  (1+cos  <j))P 
d<+)    K 

4>-»-180 

Evaluating 

lim  (sin  $)1_P  (1+cos  $)P  =  lim  (sin  <j>)  (-£§p$-) 
$-►180  <}>->-180 

by  applying  1'  Hopitale  rule  to  the  term 


.  .   1+cos  t>   _  ,  .  -sin 
lim  1 —  =  lim 


sin  1>        cos  $ 
0-180         (J)->18  0 

Thus,  the  product  sin  -M^^nr)*  approaches  zero  for 
any  p>l,  since  both  terms  in  this  product  approach  zero 
as  $+180. 
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Figure  39  shows  a  plot  of  computed  values  of  JcJk/V  | 
for  different  values  of  p>l.   From  this  we  can  draw  the 
following  summarizing  conclusions. 
For  l<p<2  the  turn  rate  has 

-a  maximum  value  at   }=cos_1 (p/2) 
-the  zero  value  at  $=0°  and  <p=180°. 
For  p=2  the  turn  rate 
-is  zero  at  <j>=180  . 

-monotonically  increasing  with  decreasing  <j> 
approaching  the  value  (4V  /K)  when  $=0° . 
For  p  2  the  turn  rate 
-is  zero  at  ^=180°. 

-is  monotonically  increasing  with  decreasing  $ 

i  •  «o 

approaching  °°  as  $  =  0  . 

As  seen  from  the  above  for  p>2,  the  turn  rate  increases 
monotonically  with  decreasing   approaching  an  infinite 
turn  rate  to  hit  a  point  target.   To  avoid  this  singularly 
we  must  make  some  provision.   If  the  torpedo's  maximum 
turn  rate  is  exceeded  only  at  some  very  small  range,  a  hit 
is  likely. 

Figure  40  illustrates  the  hit  criterion  used.   Assume 
a  rudder  of  length  L_  is  situated  directly  behind  the 
propeller  which  is  idealizer  as  an  acoustic  point  source. 

If  the  torpedo  becomes  turn  rate  limited  at  some  range 
r=r',  it  at  best  can  proceed  along  a  circular  path  which 
lags  the  desired  pursuit  trajectory,  or  at  worst  it  can 
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loose  acoustic  contact  and  go  into  "hold-in,"  maintaining 

a  constant  heading  at  the  angle  <j> '  .   In  this  last  case,  it 

will  cross  the  line  of  target  ship's  course  in  the  time: 

t  =  r'/V  (118) 

r 

If  this  advance  is  less  than  the  length  of  the  rudder  LR, 
the  torpedo  will  impact  the  rudder.   Kence ,  the  limiting 
condition  is 

Si-  L 

p     ^R 


r«  =  PLR  (119) 

For  this  analysis,  we  have  arbitrarily  chosen 

L   =  3  m  (120) 

R 

Thus,  the  torpedo's  turn  rate  has  not  been  exceeded 
when  the  range  to  the  target  is 

r'  =  pLR  =  3  p  (121) 

and  a  hit  is  assumed. 

Now  we  can  analyze  the  cases  p>2  and  Kp<2  on  the  same 
footing.   If  we  match  the  torpedo's  maximum  turn  rate  to 
a  particular  range  r',  it  is  certain  that  its  turn  rate 
has  not  been  exceeded  earlier  in  the  run. 

This  turn  rate-range  matching  is  done  by  substituting 
the  "hit-range"  defined  as  r'  =  3  p  into  Eq.  (106)  and 
solving  it  for  sin  $,  yielding 
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Vs 

W  =  -  if  sin  * 


^ax(3p) 

sin  i  =  - • • 


V 
s 


Since  J)  is  always  negative,  sin  <$>   is  always  positive  and 
equal  to 

3p|i| 

.:„  .  _  '  'max 

Sln  * v (122) 

s 

Two  values  of  $  satisfies  this  equation,  and  they  are  denoted 


3p 
PA  =   sin     — '  T'   *;  for  <J>  <  90"  (123) 


=  sin"1   -ZJllmax.  for  ^  <  9Qo 


s 


*B  =  180°  "  *A 


A  unique  value  of  K  may  now  be  found  using  Eq.  (Ill)  for 
each  of  the  angles  $   and  $    . 

A        B 


V  (1+cos  4>)P 
K  =  -= 


*Uax(sin  ^P"2 


These  values  of  K  are  designated  as  K   and  K  ,  and  for 
each  there  is  a  corresponding  value  of  <j>   from  Eq.  (109) 
with  o      given  by  the  initial  detection  range 

(*0>Aand  (*0»B 
In  the  case  of  l<p<2,  the  turn  rate  does  not  necessarily 
increase  monotonically  during  the  pursuit  homing  trajectory, 
and  we  must  check  that  either  <t>   or  4>  ,  respectively,  is  not 
beyond  the  angle  corresponding  to  3>   x» 

95 


As  seen  from  Eq.  (113) ,  the  turn  rate  reaches  a  maximum 
at  an  angle  given  by 

$*   =   cos"1(P/2) 
with  a  corresponding  turn  rate  given  by 

V  (1+P/2)P  [1-(P/2)2]1-P/2 
1  ' max  K 

Thus,  the  value  of  K  =  K*  for  which  the  limiting  turn 

rate  is  achieved 

(1+P/2)P[1-(P/2)2]1"P/2 


K  =  *- 


'  '  max 
It  is  important  to  note  that,  if  $-<<$>*<  (<J>  ),,  the  "A" 
trajectory  is  invalid  since  the  maximum  turn  rate  of  the 
torpedo,  reached  at  the  range  r  =  3p  at  4>  =4>   is  exceeded 
earlier  in  the  trajectory.   In  that  case,  the  limiting 
trajectory  is  the  "*"  trajectory.   Along  the  same  lines, 
we  argue  that  the  "B"  trajectory  always  is  a  limiting 
trajectory,  since  the  equality  <J>_«P*  (^0)B  cannot  be 
satisfied.   This  follows  from  the  fact  that 

0°<V*  <   90° 
and  ?   90°  as  ?   is  the  supplement  of  <J>A« 

Below  some  critical  ship  speed,  the  torpedo  will  not 
be  turn  rate  limited.   This  speed  is  obtained  when  the 
trajectory  is  normal  to  the  ships  velocity  vector,  ^  =90  , 
at  a  range  of  r  =  3p. 
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Thus,    from  Eq.     (106),    we   obtain 

<*>I*L«  =  <vs>nl  sin  90° 

yielding 

(Vnl  "  l3P'lilmax  (124) 

Then  introducing  p  =  V  / (V  ) MT ,  we  get 

<Vs>NL  =  ^V'^max^  <125> 

Furthermore,  we  see  from  Eq.  (108)  that  for  4>=90° 

(K)NL  =  3p  (126) 

Again,  by  probing  Eq.  (88)  we  can  obtain  the  corres- 
ponding values  of  <P    ,  designated  {<$>    )   ,  where  the  subscript 
"NL"  is  used  to  indicate  the  "no-limit"  boundary  point. 

The  computational  procedure  is  based  on: 

1.  Assuming  a  torpedo  speed  V  =  35  kts. 

2.  Using  an  initial  detection  range  (beginning  of 
homing)  r   =  250  m. 

3.  Assuming  the  following  maximum  turn  rates: 
a.   8  o/s    b.   12  o/s     c.   16  o/s 

d.   24  o/s   e.   36  o/s     f.   48  o/s 

4.  The  following  range  of  target  speeds 

0<V  <25  kts 

s 

The  aim  of  the  computation  is  to  determine  whether  the 
limiting  £  ,  and  hence  AOB,  is  governed  by  <j>*  or  by  the 
turn  rate  at  the  range  r  =  3p. 
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The  calculations  are  devided  into  two  parts,  and  are 
performed  on  a  TI  59  calculator. 

For  a  given  ship  speed,  Part  I  gives  the  sequential 
calculations  of  (J>  ,  K  ,  <$>    ,    K  ,  0*  and  K*  for  each  of  the 
turn  rates.   The  computer  program  is  given  in  Appendix  E. 

Then,  in  Part  II,  the  probe  calculation  for  (<}>)*/ 

(*o}B'  (5o}*  =  (:po}A  B  *  are  Performed-   The  program  is 
given  in  Appendix  E. 

The  limiting  results  are  given  on  a  polar  plot,  Fig.  41 
with 

(A0B)A,B,«  "  180°  "  (»0>A,B,* 

As  seen  from  Fig.  41,  for  a  target  speed  of  15  kts,  we 

need  AOB>85°  at  <b   for  a  maximum  turn  rate   16  o/s  in  order 

o 

to  have  a  hit  at  the  first  attack. 
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IX.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  detection  performance  of  a  passive  homing  torpedo 
used  against  shallow-draft  surface  ships  operating  in 
moderate  sea  states  was  investigated.   Attention  was 
focused  on  the  effects  of  scattering  from  the  randomly- 
rough  sea  surface  and  scattering  and  absorption  from  the 
bubble-dominated  inhomogeneous  layer  just  below  the  sea 
surface.   The  effects  of  these  two  scattering  mechanisms 
were  separately  estimated  and  their  relative  importance 
were  compared. 

The  passive  sonar  equation  was  used  to  predict  the 
performance  of  the  homing  system,  and  the  detection  range 
considering  these  two  scattering  effects  was  obtained  and 
compared  to  the  detection  range  based  on  a  reference 
model.   An  idealized  propagation  model  was  used  as 
reference  of  comparison.   This  reference  model  was  based 
on  a  noise  source  model  for  a  cavitating  propeller,  the 
operational  characteristics  for  a  square-law  detector, 
and  a  transmission  model  associated  with  a  homogeneous, 
unbounded  medium. 

Due  to  high  frequency,  moderate  sea  state  and  low 
grazing  angles,  the  scattering  from  the  randomly  rough 
sea  surface  was  found  to  be  dominated  by  the  direct  path. 
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This  result  includes  effects  from  geometrical 
shadov/ing. 

The  effect  of  scattering  and  absorption  from  the 
bubble-dominated,  inhomogeneous  subsurface  layer  was  in- 
vestigated using  multiple  scattering  theory.   Both  the 
coherent  and  incoherent  limits  were  investigated  by  incor- 
porating the  associated  absorption  coefficient  into  the 
transmission  equation. 

The  effect  on  the  sound  speed  from  the  bubble  content 
was  found  to  be  negligible. 

At  the  assumed  depth  setting  of  6  m  for  the  torpedo's 
search  and  attack  phase,  the  scattering  from  the  bubbles 
increased  the  transmission  loss.   This  increase  depended 
on  the  frequency  and  the  wind  speed. 

Two  operating  frequencies  were  investigated,  60  kHz 
and  30  kHz.   For  both  cases,  bubbles  significantly 
decreased  the  detection  range. 

For  a  torpedo  system  operating  at  the  high  frequencies, 
e.g.,  60  kHz,  the  result  indicates  the  inadvisability  of 
using  a  searching  and  attack  depth  near  the  surface,  i.e., 
z  =  6  m. 

For  an  operating  frequency  around  30  kHz,  the  calculated 
detection  ranges  is  such  that  the  turn  rate  requirements 
for  a  pursuit  homing  trajectory  become  excessive.   For  a 
maximum  turn  rate  of  16  o/s,  this  limitation  can  be  avoided 
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by  adapting  a  tactical  procedure  where  the  angle  on  bow 
(AOB)  at  the  beginning  of  the  torpedo  attack  is  greater 
than  85°. 

At  sea  state  3,  the  results  show  a  consistent  and 
general  trend  towards  the  need  for  lower  operating  fre- 
quency in  order  to  increase  the  detection  range.   An 
operating  frequency  below  30  kHz  seems  indicated. 

Furthermore,  a  search  depth  below  the  bubble-dominated 
subsurface  layer,  i.e.,  z>15-20  m  would  result  in  an 
increased  detection  range. 

To  reduce  the  operational  limitations  induced  by  the 
scattering  and  absorption  effects,  a  high  maximum  turn 
rate  together  with  a  variable  speed  capability,  where 
p-2  would  be  beneficial. 

The  result  of  this  analysis  has  clearly  demonstrated 
the  importance  of  environmental  factors  on  the  torpedo 
capability,  and  is  useful  in  giving  insight  into  the  be- 
havior of  a  homing  torpedo  during  its  critical  attack 
phase. 

A  valuable  follow-on  of  this  study  would  be  an  investi- 
gation of  the  effects  of  the  bubble-dominated  subsurface 
layer  on  target  validation  and  pitch  plane  steering  when 
the  torpedo  search  depth  is  50-60  m. 
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APPENDIX  A 
DETAILED  OCEANOGRAPHIC  BACKGROUND  MATERIAL 

A.  GEOGRAPHIC  DESCRIPTION 

The  Norwegian  coastal  waters  constitute  the  eastern 
boundary  of  the  Norwegian  Sea.   Although  some  general 
aspects  related  to  the  Norwegian  Sea  will  be  covered, 
this  analysis  will  be  concentrated  on  the  Norwegian 
coastal  waters  above  68°N. 

B.  WIND 

The  northern  region  of  the  Norwegian  Sea  is  affected 
by  the  Polar  Easterlies  and  the  southern  region  by  the 
prevailing  Westerlies.   There  are  two  dominant  air  masses 
which  are  relatively  permanent: 

-The  Greenland  high. 

-The  Iceland  low. 
These  pressure  systems  produce  storms  which  are  carried 
across  the  Norwegian  Sea  in  a  belt  from  Iceland  towards 
the  Norwegian  Coast  causing  steady  precipitation  and  wind 
most  of  the  year.   The  steep  Norwegian  Coast  has  a  con- 
siderable influence  on  the  winds  and  consequently  also  on 
the  waves  in  the  coastal  waters.   The  main  general  modi- 
fications are  that  the  streamlines  tend  to  run  parallel 
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to  the  coast  and  that  wind  and  sea  increases  with  distance 
from  the  sheltered  coast  into  open  ocean. 

Strong  local  variation  may  occur.   The  most  important 
of  these  are  the  marked  local  increase  in  wind  speed  in 
areas  where  the  coast  sharply  changes  direction.   One  such 
"corner  effect,"  caused  by  the  confluence  of  the  streamlines, 
occurs  near  "Nordkapp"  (North  Cape) .   Also  of  importance 
are  monsoonal  effects  due  to  the  different  heat  capacities 
between  the  ocean  and  the  continent.   Drainage  of  cold  air 
from  the  inland  valleys  in  the  wintertime  causes  a  marked 
increase  in  the  wind  speed  in  several  areas  along  the 
coast.   Most  of  these  coastal  effects  are  significantly 
dissipated  at  distances  of  approximately  50  nmi.  from  the 
coast.   The  fact  that  the  wind  tends  to  blow  along  the 
coast  is  clearly  demonstrated  in  Figs.  12  and  13  which 
include  data  from  weather  stations  from  "Hillesoy"  to 
"Ona"  and  "Myken"  to  "Furuholmen"  respectively.   The  high 
frequency  of  offshore  winds  is  caused  by  the  drainage  of 
cold  air  from  the  inland  valleys  during  winter  time. 
That  this  phenomenon  is  closely  connected  to  the  coast 
is  illustrated  by  the  fact  that  it  is  missing  at  weather 
station  "Skomvaer"  situated  approximately  50  nmi.  off  the 
main  coast.   A  frequency  distribution  of  observed  wind 
speeds  along  the  coast,  obtained  from  Ref.  1,  is  presented 
in  Table  V.   A  summary  of  this  table  follows: 
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TABLE  V 

FREQUENCY  DISTRIBUTION  OF  WIND  SPEED  IN  PERCENT  PER  YEAR 
AT  WEATHER  STATIONS  ALONG  THE  NORWEGIAN  COAST 


m/sec 

0-7 

8-13 

14-20 

$2.1 

Beaufort 

0-4 

5+6 

7+8 

^9 

Ferder 

67 

28 

5 

0.1 

Lyngor 

79 

19 

2 

0.2 

Lista 

67 

28 

5 

0.2 

Utsira 

73 

22 

5 

0.3 

Hellisoy 

77 

19 

4 

0.2 

Krakenes 
Ona 

58 

70 

28 
22 

11 
7 

2.7  Jan-Dec 
0.7  1^49-1975 

Sula 

61 

30 

8 

1.1 

Nordoyan 

49 

35 

14 

1.9 

Myken 

64 

26 

9 

1.2 

Skomvaer 

58 

31 

10 

0.9 

Andenes 

79 

18 

3 

0.1 

Torsvar 

72 

23 

5 

0.4 

Fruholmen 

54 

32 

12 

1.8 

Vardo 

75 

22 

3 

0.1 

Biornoya 

63 

31 

6 

0.4  1956-1975 

Polarfront 

46 

30 

14 

1.4 
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The  highest  winds  are  reported  from  the  areas  between 
62  N  and  68  N.   The  frequency  of  high  winds  in  this  area 
are  significantly  greater  than  those  found  at  stations  to 
the  south  of  62  N  and  also  to  the  north  of  68°N.   The  wind 
condition  around  North  Cape,  are  very  severe.   This  area 
can  be  compared  with  the  other  coastal  area  of  high  wind 
speed  such  as  "Stadt"  (represented  by  the  weather  station, 
"Krakenes" ) .   In  these  areas  the  frequency  of  storms  is 
greater  than  at  "Polarfront"  situated  in  the  open  Norwegian 
Sea. 

C .   WAVES 

Frequency  distribution  of  significant  wave  heights  are 
represented  in  Table  VI. 

The  station  north  of  68  N  is  characterized  by  compara- 
tively small  frequency  of  high  waves.   Even  at  "Furuholmen" 
where  wind  conditions  are  very  severe,  the  frequency  dis- 
tribution of  significant  wave  heights  is  similar  to  more 
sheltered  areas  like  "Utsira." 

The  seasonal  variation  for  the  area  of  interest  around 
70°N  is  given  in  Figs.  3  and  4.   The  average  monthly  dis- 
tribution of  significant  wave  heights  for  a  typical 
station  like  "Andenes"  is  given  in  Fig.  3   and  Fig.  4 
gives  the  yearly  distribution  for  this  station. 
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TABLE  VI 

FREQUENCY  DISTRIBUTION  OF  SIGNIFICANT  WAVE 
HEIGHT  IN  PERCENT  AT  THE  WEATHER  STATION  "ANDENES" 


STATION:   Andenes  YEARS:  1949-1972 

SS     0+1+2  3  4  5  6  7  8  9 

HO-  0.5  -  1.3  -  2.5-4.0  6.0  -  9.0-  13.9   N 

0.4  1.2  2.4  3.9  5.9  8.9  13.9 

JAN    15.3  26.9  33.3  15.5  6.6  2.3  0.18  0.00   2542 

FEB    17.7  28.7  31.3  14.2  5.9  2.0  0.14  0.00   2295 

MAR    18.1  28.9  30.8  14.3  5.8  1.9  0.17  0.00   2529 

APR    24.5  33.1  27.4  10.8  3.4  0.8  0.04  2293 

MAY    30.9  35.5  23.9  7.7  1.7  0.3  0.01  1720 

JUN    37.4  36.9  20.2  4.8  0.6  0.1  0.00  1565 

JUL    42.6  36.2  17.1  3.7  0.4  0.0  0.00  1618 

AUG    40.6  36.4  18.0  4.3  0.5  0.1  1621 

SEP    29.7  34.3  24.4  8.9  2.1  0.5  0.01  2301 

OCT    21.8  31.3  29.4  12.4  4.1  1.0  0.01  2374 

NOV    21.7  31.5  29.9  11.7  4.1  1.1  0.06  2303 

DEC    18.5  30.1  32.2  13.2  4.7  1.2  0.09  2379 

YEAR   26.6  32.5  26.5  10.1  3.3  0.9  0.06  0.00 


N   =  Number  of  Observations 

SS   =  Class  Interval,  State  of  Sea 

Hq  =  Significant  Wave  Height 
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D.   AMBIENT  NOISE 

Due  to  the  fact  that  the  Norwegian  Sea  is  physically 
separated  from  the  Atlantic  by  the  Faeroy-Shetland-Iceland 
ridge,  and  from  the  Greenland  Sea  by  the  Jan  Mayan  ridge, 
little  long  distance  shipping  noise  is  transferred  into 
the  area.   This,  combined  with  relatively  low  shipping 
traffic  in  the  central  and  northern  parts  of  the  Norwegian 
Sea,  produces  a  relatively  low  ambient  noise  level  for  the 
frequency  band  100<f<1000  Hz,  especially  when  noise  from 
marine  life  is  not  included.   Further,  in  the  absence  of 
nearby  shipping  and  marine  life,  the  ambient  noise  level 
in  the  frequency  band  Kf<50  kHz  is,  according  to  Ref.  4, 
dominated  by  the  wind. 

A  typical  area  for  the  central  part  of  the  Norwegian 
Sea  can  be  represented  by  the  weather  station  "Polarf ront" 
at  66  N,  2  E.   Reference  1  shows  that  there  is  approximately 
a  15%  chance  of  finding  wind  forces  of  Beaufort>6.   These 
effects  predict  a  moderate  ambient  noise  level  in  the 
frequency  range  100  Hz-50  kHz  for  the  central  and  northern 
regions  of  the  Norwegian  Sea. 

The  shallow  coastal  Norwegian  waters  are,  according  to 
Ref.  4,  typically  5-10  dB  noisier  than  the  corresponding 
deep  waters.   However,  great  variability  caused  by  local 
traffic,  fishing  fleet  activity,  marine  life,  and  local 
wind  conditions  makes  ambient  noise  level  prediction 
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difficult  in  these  areas.   This  means  that  accurate  ambient 
noise  level  determinations  have  to  be  made  on  the  spot,  as 
it  is  both  site  and  time  dependent. 

A  noticeable  influence  on  the  ambient  noise  level  is 
rain,  which  is  a  year  around  feature  along  the  Norwegian 
coast.   As  seen  from  Fig.  8,  taken  from  Ref.  4,  rain  has 
a  tendency  to  produce  a  constant  high  ambient  noise  level 
over  a  large  frequency  range,  thus  dominating  other 
effects.   Furthermore,  for  the  upper  frequency  of  interest, 
i.e. ,  around  60  kHz  the  lower  bound  for  the  ambient 
noise  is  determined  by  the  thermal  agitation,  see  Fig.  8. 

In  determining  the  figure  of  merit  (FOM)  for  a 
passive  sonar  system,  the  noise  level  will  be  the  larger 
of  either  the  self  noise  or  the  ambient  noise.   For  a 
torpedo  the  self  noise  will  typically  be  dominant. 

E.   SOUND  SPEED  PROFILES 

Again,  concentrate  on  data  relevant  to  Norwegian 
coastal  waters.  According  to  Ref.  6,  which  covers  the 
southern  part  of  the  Norwegian  coast,  low  sound  speeds 
are  common  because  of  the  influence  of  water  from  the 
Baltic  Sea  combined  with  fresh  water  drainage  from  the 
fjords.  Furthermore,  great  variability,  both  seasonal 
and  within  seasons,  is  encountered.  Figure  9  obtained 
from  Ref.  5  gives  a  picture  of  the  sound  speed  profiles 
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for  the  northern  Norwegian  coast.   Again,  large  variations 
are  common.   Noticeable  in  both  sets  of  data  is  a  typical 
seasonal  pattern  of  strong  cooling  of  the  surface  layer 
during  winter  and  a  similarly  strong  heating  during 
summer .   Furthermore ,  note  that  the  minimum  and  maximum 
are  relatively  shallow,  i.e.,  less  than  50  m. 

Also  characteristic  is  the  influence  of  the  cold  and 
fresh  melt  waters  drained  out  through  the  fjord-arms 
during  spring  and  summer. 

To  illustrate  the  sonar  problems  associated  with  these 
sound  speed  profiles,  ray  paths  for  the  extremes  of  Fig.  9 
are  shown  in  Figs.  10  and  11,  where  the  source  is  3  m  below 
the  sea  surface. 
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APPENDIX  B 
SURFACE  SCATTERING  TI  59  PROGRAM 

A.   INTRODUCTION 

This  program  gives  specular  scattered  power  at  the 
receiver  versus  incoming  power  at  the  randomly  rough 
surface  in  the  high  frequency  limit  according  to  Eq.  (70) , 

Shadowing  of  surface  areas  by  other  parts  of  the 
boundary  are  taken  into  account  by  the  bistatic  shadowing 
function  S ( 9) . 

Furthermore,  the  program  gives  the  effect  of  the 
randomly  rough  surface  compared  to  the  idealized  free- 
field  condition  as  expressed  in  Eq.  (80) .   The  results 
of  these  calculations  is  given  in  Fig.  27. 

3.   PROGRAM  STEPS 

A  block-diagram  of  the  computer  program  is  given  in 
Fig.  26. 

The  program  uses  the  partitioning  ratio  of  program 
to  data  space  according  to  code  4  OP17.  The  users  in- 
structions are  as  follows: 
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Procedure 


Enter 


Press 


Display 


Enter  data 

Enter  data 

Enter  data 

Enter  data 

Enter  data 

Calculate  9, 
Ri ,Rz,A 

Calculate  <;12> 

Calculate  v 

Calculate  erfc  v 

Calculate  S(9) 


<s2> 


Detection  range 
Source  depth 
Receiver  depth 
Beam  width 
Wind  speed 


Calculate  <p1pl*> 
Calculate  AIL 


2nd  A 

V 

R/S 

h  l 

s 

R/S 

h  1 
r 

R/S 

t1 

R/S 

w1 

2nd  B 

A 

(9,Ri,R2,A) 

2nd  C 

<?12>1 

2nd  D 

v1 

A 

.    1 
erfc  v 

B 

SC9)1 

<s2>1 

C 

<PiPi*> 

D 

AIL1 

'These  values  are  printed  automatically  if  the  calculator 
is  connected  to  the  PC-100A  Print  Cradle. 


For  the  error  function  complement  we  have 


er 


fc(v)  =   —     f    e"a  da 
/¥  v 


in  contrast  to  the  normal  distribution 


Q(u)  =  —  /  e 

/2rf  a 


-t2/2 


dt 


However,  there  is  a  linear  relationship  between  the 

two  functions. 

Ill 


The  numerical  equation  used  to  calculate  the  erfc(v) 
is  a  modified  program  from  Texas  Instruments  [Ref.  24]. 
erfc(v)  =  Z  (v)  [bia+bac^+bac^+b^+bsa5] 


where 


a,   = 


1+pv 

p   =  .231649 

b   =  .451673691 
i 

b2  =  -.504257336 
b3  =  2.51939026 
b„  =  -2.563346623 
b5  =  1.881292139 

The  program  steps  are  listed  below,  giving  location 
(LOC) ,  code  (COD) ,  key  symbol  (KEY) ,  and  comments. 


000 

76 

LBL 

001 

16 

A' 

002 

42 

ST 

003 

00 

00 

004 

99 

PRT 

005 

91 

R/S 

006 

42 

STO 

007 

01 

01 

008 

99 

PRT 

009 

91 

R/S 

010 

42 

STO 

011 

02 

02 

012 

99 

PRT 

013 

91 

R/S 

014 

42 

STO 

015 

03 

03 

016 

99 

PRT 

017 

91 

R/S 

018 

42 

STO 

019 

04 

04 

020 

99 

PRT 

021 

91 

R/S 

022 

76 

LBL 

023 

17 

B« 

024 

53 

( 

025 

53 

( 

026 

43 

RCL 

027 

00 

00 

028 

65 

X 

029 

43 

RCL 

030 

02 

02 

031 

54 

) 

032 

55 

* 

033 

53 

( 

034 

43 

RCL 

035 

01 

01 

036 

85 

+ 

037 

43 

RCL 

038 

02 

02 

039 

54 

) 

040 

54 

) 

041 

42 

STO 

042 

05 

05 

043 

53 

( 

044 

53 

( 

045 

43 

RCL 

046 

05 

05 

047 

55 

T 

048 

43 

RCL 

049 

02 

02 

050 

54 

) 

051 

22 

INV 

052 

30 

TAN 

053 

54 

) 

054 

42 

STO 

055 

06 

06 

056 

99 

PRT 

057 

53 

( 

058 

53 

( 

059 

43 

RCL 

060 

05 

05 

061 

55 

V 

062 

53 

( 

063 

43 

RCL 

064 

06 

06 

065 

38 

SIN 

066 

54 

) 

067 

54 

) 

068 

42 

STO 

112 


069 

07 

07 

070 

99 

PRT 

071 

53 

( 

072 

53 

( 

073 

43 

RCL 

074 

00 

00 

075 

75 

- 

076 

43 

RCL 

077 

05 

05 

078 

54 

) 

079 

55 

T 

080 

53 

( 

081 

43 

RCL 

082 

06 

06 

083 

38 

SIN 

084 

54 

) 

085 

54 

) 

086 

42 

STO 

087 

08 

08 

088 

99 

PRT 

089 

53 

( 

090 

43 

RCL 

091 

07 

07 

092 

85 

+ 

093 

43 

RCL 

094 

08 

08 

095 

54 

) 

096 

42 

STO 

097 

09 

09 

098 

25 

CLR 

099 

69 

OP 

100 

00 

00 

101 

01 

1 

102 

03 

3 

103 

03 

3 

104 

05 

5 

105 

01 

1 

106 

07 

7 

107 

01 

1 

108 

03 

3 

109 

69 

OP 

110 

04 

04 

111 

53 

( 

112 

53 

( 

113 

89 

TT 

114 

65 

X 

115 

43 

RCL 

116 

08 

08 

117 

54 

) 

118 

65 

X 

119 

53 

( 

120 

43 

RCL 

121 

03 

03 

122 

38 

SIN 

123 

65 

X 

124 

43 

RCL 

125 

07 

07 

126 

54 

) 

127 

54 

) 

128 

42 

STO 

129 

10 

10 

130 

69 

OP 

131 

06 

06 

132 

91 

R/S 

133 

76 

LBL 

134 

18 

c 

135 

25 

CLR 

136 

69 

OP 

137 

00 

00 

138 

00 

00 

139 

06 

6 

140 

02 

2 

141 

07 

7 

142 

03 

3 

143 

02 

2 

144 

03 

3 

145 

03 

3 

146 

69 

OP 

147 

04 

04 

148 

53 

( 

149 

53 

( 

150 

53 

( 

151 

43 

RCL 

152 

04 

04 

153 

65 

X 

154 

05 

5 

155 

93 

• 

156 

01 

1 

157 

02 

2 

158 

54 

) 

159 

85 

+ 

160 

03 

3 

161 

54 

) 

162 

65 

X 

163 

93 

• 

164 

00 

0 

165 

00 

0 

166 

01 

1 

167 

54 

) 

168 

50 

lxl 

169 

42 

STO 

170 

11 

11 

171 

69 

OP 

172 

06 

06 

173 

91 

R/S 

174 

76 

LBL 

175 

19 

D' 

176 

53 

( 

177 

43 

RCL 

178 

06 

06 

179 

30 

TAN 

180 

35 

1/x 

181 

54 

) 

182 

42 

STO 

183 

12 

12 

184 

53 

( 

185 

43 

RCL 

186 

12 

12 

187 

55 

T 

188 

53 

( 

189 

53 

( 

190 

02 

2 

191 

65 

X 

192 

43 

RCL 

193 

11 

11 

194 

54 

) 

195 

34 

& 

196 

54 

) 

197 

54 

) 

198 

42 

STO 

199 

13 

13 

200 

99 

PRT 

201 

91 

R/S 

202 

76 

LBL 

203 

11 

A 

204 

53 

( 

205 

53 

( 

206 

53 

( 

207 

53 

( 

208 

43 

RCL 

209 

13 

13 

210 

33 

X2 

211 

54 

( 

212 

55 

f 

213 

02 

2 

214 

54 

) 

215 

94 

+/- 

216 

22 

INV 

217 

23 

LNX 

218 

54 

) 

219 

65 

X 

220 

53 

( 

221 

53 

( 

113 


222 

89 

77 

223 

34 

/x 

224 

65 

X 

225 

02 

2 

226 

54 

) 

227 

35 

l/x 

228 

54 

) 

229 

54 

) 

230 

42 

STO 

231 

14 

14 

232 

53 

( 

233 

53 

( 

234 

53 

( 

235 

53 

( 

236 

43 

RCL 

237 

13 

13 

238 

65 

X 

239 

93 

. 

240 

02 

2 

241 

03 

3 

242 

01 

1 

243 

06 

6 

244 

04 

4 

245 

01 

1 

246 

09 

9 

247 

86 

+ 

248 

01 

1 

249 

54 

) 

250 

35 

l/x 

251 

54 

) 

252 

42 

STO 

253 

15 

15 

254 

53 

( 

255 

53 

( 

256 

43 

RCL 

257 

14 

14 

258 

65 

X 

259 

01 

1 

260 

93 

. 

261 

08 

8 

262 

01 

1 

263 

54 

) 

264 

65 

X 

265 

53 

( 

266 

43 

RCL 

267 

15 

15 

268 

65 

X 

269 

93 

. 

270 

03 

3 

271 

01 

1 

272 

09 

9 

273 

03 

3 

274 

08 

8 

275 

01 

1 

276 

05 

5 

277 

03 

3 

278 

75 

- 

279 

93 

. 

280 

03 

3 

281 

05 

5 

282 

06 

6 

283 

05 

5 

284 

06 

6 

285 

03 

3 

286 

07 

7 

287 

08 

8 

288 

02 

2 

289 

65 

X 

290 

43 

RCL 

291 

15 

15 

292 

33 

x2 

293 

85 

+ 

294 

01 

1 

295 

93 

. 

296 

97 

7 

297 

08 

8 

298 

01 

1 

299 

04 

4 

300 

07 

7 

301 

07 

7 

302 

09 

9 

303 

03 

3 

304 

07 

7 

305 

65 

X 

306 

43 

RCL 

307 

15 

15 

308 

45 

Yx 

309 

03 

3 

310 

75 

- 

311 

01 

1 

312 

93 

• 

313 

08 

8 

314 

02 

2 

315 

01 

1 

316 

02 

2 

317 

05 

5 

318 

05 

5 

319 

09 

9 

320 

07 

7 

321 

08 

8 

322 

00 

0 

323 

65 

X 

324 

43 

RCL 

325 

15 

15 

326 

45 

Yx 

327 

04 

4 

328 

85 

+ 

329 

01 

1 

330 

93 

. 

331 

03 

3 

332 

03 

3 

333 

00 

0 

334 

02 

2 

335 

07 

7 

336 

04 

4 

337 

04 

4 

338 

02 

2 

339 

09 

9 

340 

00 

0 

341 

65 

X 

342 

43 

RCL 

343 

15 

15 

344 

45 

Yx 

345 

05 

5 

346 

54 

) 

347 

54 

) 

348 

42 

STO 

349 

16 

16 

350 

99 

PRT 

351 

91 

R/S 

352 

76 

LBL 

353 

12 

B 

354 

53 

( 

355 

53 

( 

356 

53 

( 

357 

43 

RCL 

358 

13 

13 

359 

33 

X2 

360 

94 

+/- 

361 

22 

INV 

362 

23 

LNX 

363 

54 

) 

364 

75 

- 

365 

53 

( 

366 

89 

TT 

367 

34 

/x 

368 

65 

X 

369 

43 

RCL 

370 

13 

13 

371 

65 

X 

372 

43 

RCL 

373 

16 

16 

374 

54 

) 
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375 

54  ) 

426 

99  PRT 

376 

55  t 

427 

91  R/S 

377 

53  ( 

428 

76  LBL 

378 

89  it 

429 

13  C 

379 

34  /x 

430 

25  CLR 

380 

65  x 

431 

69  OP 

381 

04  4 

432 

00  00 

382 

65  x 

433 

03  3 

383 

4  3  RCL 

434 

06  6 

384 

13  13 

435 

01  1 

385 

54  ) 

436 

05  5 

386 

54  ) 

437 

01  1 

387 

42  STO 

438 

03  3 

388 

17  17 

439 

03  3 

389 

53  ( 

440 

07  7 

390 

53  ( 

441 

69  OP 

391 

53  ( 

442 

04  04 

392 

01  1 

443 

53  ( 

393 

75  - 

444 

53  ( 

394 

53  ( 

445 

43  RCL 

395 

53  ( 

446 

10  10 

396 

53  ( 

447 

65  x 

397 

4  3  RCL 

448 

43  RCL 

398 

17  17 

449 

18  18 

399 

65  x 

450 

54  ) 

400 

04  4 

451 

55  t 

401 

54  ) 

452 

53  ( 

402 

94  +/- 

453 

53  ( 

403 

54  ) 

454 

53  ( 

404 

22  INV 

455 

53  ( 

405 

23  LNX 

456 

4  3  RCL 

406 

54  ) 

457 

07  07 

407 

54  ) 

458 

33  x2 

408 

65  x 

459 

54  ) 

409 

53  ( 

460 

65  x 

410 

01  1 

461 

08  8 

411 

75  - 

462 

54  ) 

412 

43  RCL 

463 

65  x 

413 

16  16 

464 

89  7T 

414 

54  ) 

465 

54  ) 

415 

54  ) 

466 

65  x 

416 

55  i 

467 

4  3  RCL 

417 

53  ( 

468 

11  11 

418 

4  3  RCL 

469 

54  ) 

419 

17  17 

470 

54  ) 

420 

65  x 

471 

42  STO 

421 

04  4 

472 

19  19 

422 

54  ) 

473 

69  OP 

423 

54  ) 

474 

06  06 

424 

42  STO 

475 

91  R/S 

425 

18  18 

476 

76  LBL 

477 

14 

D 

478 

25 

CLR 

479 

69 

OP 

480 

00 

00 

481 

01 

1 

482 

06 

6 

483 

01 

1 

484 

04 

4 

485 

69 

OP 

486 

04 

04 

487 

53 

( 

488 

53 

( 

489 

43 

RCL 

490 

19 

19 

491 

28 

LOG 

492 

54 

) 

493 

65 

X 

494 

01 

1 

495 

00 

0 

496 

54 

) 

497 

42 

STO 

498 

20 

20 

499 

69 

OP 

500 

06 

06 

501 

91 

R/S 
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APPENDIX  C 
BUBBLE  DYNAMICS  TI  59  PROGRAM 

A.   INTRODUCTION 

This  program  gives  the  resonant  frequency  f   and  the 
damping  constant  5  for  bubbles  according  to  Eq.  (92) . 
Furthermore,  the  program  gives  scattering  cross  section 
a  s,  extinction  cross  section  a  ,  and  absorption  cross 
section  a^  as  a  function  of  bubble  radius  a,  incoming  fre- 
quency f,  and  depth  z.   The  results  of  the  calculations 
are  given  in  Figs.  31  and  32. 

3.   PROGRAM  STEPS 


The  user's  instructions  are  as  follows: 


Procedure 


Enter 


Press 


Display 


Enter  data 

Enter  data 

Enter  data 

Calculate  X 

Calculate 
coshX  and 
sinhX 


Bubble  Radius 
Incoming  Frequency 
Depth 


Calculate 

d/b 

— 

Calculate 

b 

— 

Calculate 

6 

-■ 

Calculate 

f 

r 

— 

Calculate 

5 

-■ 

Calculate 

a 

— 

Calculate 

a 

and  a 
a 

s 

2nd  A 

a 

R/S 

f1 

R/S 

z1 

2nd  B 

X1 

2nd  C 

coshX 
sinhX1 

2nd  D 

d/b1 

2nd  E 

b1 

A 

s1 

B 

f  1 
r 

C 

61 

D 

a   1 

Oa1 
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These  values  are  printed  automatically  if  the  calculator 

is  connected  to  PC-100A  Print  Cradle. 

The  program  steps  are  listed  below  giving  location  (LOC) , 

code  (COD) ,  key  (KEY)  and  comments. 

000  76  LBL  043  55  t  086  01  01 

001  16  A'  044  43  RCL  087  99  PRT 

002  42  STO  045  11  11  088  91  R/S 

003  11  11  046  54  )  089  76  LBL 

004  99  PRT  047  85  +  090  18  C 

005  91  R/S  048  93  .  091  53  ( 

006  42  STO  049  00  0  092  53  ( 

007  12  12  050  00  0  093  43  RCL 

008  99  PRT  051  01  1  094  01  01 

009  91  R/S  052  02  2  095  22  INV 

010  42  STO  053  09  9  096  23  LNX 

011  13  13  054  65  x  097  85  + 

012  99  PRT  055  53  (  098  43  RCL 

013  91  R/S  056  01  1  099  01  01 

014  76  LBL  057  85  +  100  94  X/- 

015  17  B'  058  93  •  101  22  INV 

016  53  (  059  01  1  102  23  LNX 

017  53  (  060  65  x  103  54  ) 

018  53  (  061  43  RCL  104  55  r 

019  53  (  062  13  13  105  02  2 

020  53  (  063  54  )  106  54  ) 

021  53  (  064  54  )  107  42  STO 

022  04  4  065  54  )  108  02  02 

023  65  x  066  65  x  109  99  PRT 

024  89  :r  067  93  110  53  ( 

025  54  )  068  02  2  111  53  ( 

026  65  x  069  04  4  112  43  RCL 

027  43  RCL  070  54  )  113  01  01 

028  12  12  071  55  i  114  22  INV 

029  54  )  072  93  •  115  23  LNX 

030  65  x  073  00  0  116  75  - 

031  53  (  074  00  0  117  43  RCL 

032  53  (  075  00  0  118  01  01 

033  93  .  076  00  0  119  94  +/- 

034  00  0  077  05  5  120  22  INV 

035  00  0  078  06  6  121  23  LNX 

036  00  0  079  54  )  122  54  ) 

037  00  0  080  34  /x  123  55  * 

038  00  0  081  65  x  124  02  2 

039  00  0  082  43  RCL  125  54  ) 

040  01  0  083  11  11  126  42  STO 

041  09  9  084  54  )  127  03  03 

042  01  1  085  42  STO  128  99  PRT 
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129 

91 

R/S 

130 

76 

LBL 

131 

17 

D' 

132 

53 

( 

133 

53 

( 

134 

53 

( 

135 

53 

( 

136 

43 

RCL 

137 

01 

01 

138 

65 

X 

139 

53 

( 

140 

43 

RCL 

141 

03 

03 

142 

85 

+ 

143 

43 

RCL 

144 

01 

01 

145 

38 

SIN 

146 

54 

) 

147 

54 

) 

148 

75 

- 

149 

02 

2 

150 

65 

X 

151 

53 

( 

152 

43 

RCL 

153 

02 

02 

154 

75 

- 

155 

43 

RCL 

156 

01 

01 

157 

39 

COS 

158 

54 

) 

159 

54 

) 

160 

55 

T 

161 

53 

( 

162 

53 

( 

163 

53 

( 

164 

43 

RCL 

165 

01 

01 

166 

33 

X2 

167 

54 

) 

168 

65 

X 

169 

53 

( 

170 

43 

RCL 

171 

02 

02 

172 

75 

- 

173 

43 

RCL 

174 

01 

01 

175 

39 

COS 

176 

54 

) 

177 

54 

) 

178 

85 

+ 

179 

53 

( 

180 

53 

( 

181 

01 

1 

182 

93 

. 

183 

02 

2 

184 

65 

X 

185 

43 

RCL 

186 

01 

01 

187 

54 

) 

188 

65 

X 

189 

53 

( 

190 

43 

RCL 

191 

03 

03 

192 

75 

- 

193 

43 

RCL 

194 

01  01 

195 

38 

SIN 

196 

54 

) 

197 

54 

) 

198 

54 

) 

199 

54 

) 

200 

65 

X 

201 

01 

1 

202 

93 

. 

203 

02 

2 

204 

54 

) 

205 

42 

STO 

206 

04 

04 

207 

99 

PRT 

208 

91 

R/S 

209 

76 

LBL 

210 

10 

E 

211 

53 

( 

212 

53 

( 

213 

53 

( 

214 

53 

( 

215 

53 

( 

216 

03 

3 

217 

93 

. 

218 

02 

2 

219 

55 

220 

43 

RCL 

221 

01 

01 

222 

54 

) 

223 

65 

X 

224 

53 

( 

225 

43 

RCL 

226 

03 

03 

227 

75 

- 

228 

43 

RCL 

229 

01 

01 

230 

38 

SIN 

231 

54 

) 

232 

54 

) 

233 

55 

T 

234 

53 

( 

235 

43 

RCL 

236 

02 

02 

237 

75 

- 

238 

43 

RCL 

239 

01 

01 

240 

39 

COS 

241 

54 

) 

242 

54 

) 

243 

85 

+ 

244 

01 

1 

245 

54 

) 

246 

35 

1/x 

247 

65 

X 

248 

53 

( 

249 

54 

( 

250 

43 

RCL 

251 

04 

04 

252 

33 

X2 

253 

54 

) 

254 

85 

+ 

255 

01 

1 

256 

54 

) 

257 

35 

1/x 

258 

54 

) 

259 

42 

STO 

260 

05 

05 

261 

99 

PRT 

262 

91 

R/S 

263 

76 

LBL 

264 

11 

A 

265 

53 

( 

266 

53 

( 

267 

53 

( 

268 

01 

1 

269 

05 

5 

270 

00 

0 

271 

55 

T 

272 

53 

( 

273 

01 

1 

274 

00 

0 

275 

01 

1 

276 

03 

3 

277 

00 

0 

278 

00 

0 

279 

00 

0 

280 

65 

X 

281 

43 

RCL 
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282 

11 

11 

283 

65 

X 

284 

53 

( 

285 

01 

1 

286 

85 

+ 

287 

93 

. 

288 

01 

1 

289 

65 

X 

290 

43 

RCL 

291 

13 

13 

292 

54 

) 

293 

54 

) 

294 

54 

) 

295 

65 

X 

296 

53 

) 

297 

01 

1 

298 

75 

- 

299 

53 

( 

300 

53 

( 

301 

04 

4 

302 

93 

. 

303 

02 

2 

304 

65 

X 

305 

43 

RCL 

306 

05 

05 

307 

54 

) 

308 

35 

1/x 

309 

54 

) 

310 

54 

) 

311 

54 

) 

312 

85 

+ 

313 

01 

1 

314 

54 

) 

315 

42 

STO 

316 

06 

06 

317 

99 

PRT 

318 

91 

R/S 

319 

76 

LBL 

320 

12 

B 

321 

53 

( 

322 

53 

( 

323 

53 

( 

324 

53 

( 

325 

04 

4 

326 

93 

327 

02 

2 

328 

65 

X 

329 

43 

RCL 

330 

05 

05 

331 

65 

X 

332 

43 

RCL 

333 

06 

06 

334 

65 

X 

335 

01 

1 

336 

00 

0 

337 

01 

1 

338 

03 

3 

339 

00 

0 

340 

00 

0 

341 

00 

0 

342 

54 

) 

343 

65 

X 

344 

53 

( 

345 

01 

1 

346 

85 

+ 

347 

93 

• 

348 

01 

1 

349 

65 

X 

350 

43 

RCL 

351 

13 

13 

352 

52 

) 

353 

54 

) 

354 

55 

v 

355 

01 

1 

356 

93 

. 

357 

00 

0 

358 

02 

2 

359 

06 

6 

360 

54 

) 

361 

50 

lxl 

362 

34 

•  X 

363 

65 

X 

364 

53 

( 

365 

53 

( 

366 

89 

TT 

367 

65 

X 

368 

02 

2 

369 

65 

X 

370 

43 

RCL 

371 

11 

11 

372 

54 

) 

373 

35 

1/x 

374 

54 

) 

375 

54 

) 

376 

42 

STO 

377 

07 

07 

378 

99 

PRT 

379 

91 

R/S 

380 

76 

LBL 

381 

13 

C 

382 

53 

( 

383 

53 

( 

384 

53 

( 

385 

93 

. 

386 

00 

0 

387 

00 

0 

388 

00 

0 

389 

00 

0 

390 

04 

4 

391 

01 

1 

392 

08 

8 

393 

08 

8 

394 

07 

7 

395 

09 

9 

396 

65 

X 

397 

43 

RCL 

398 

12 

12 

399 

65 

X 

400 

43 

RCL 

401 

11 

11 

402 

54 

) 

403 

85 

+ 

404 

53 

( 

405 

43 

RCL 

406 

04 

04 

407 

65 

X 

408 

53 

( 

409 

53 

( 

410 

43 

RCL 

411 

07 

07 

412 

55 

T 

413 

43 

RCL 

414 

12 

12 

415 

54 

) 

416 

33 

X2 

417 

54 

) 

418 

54 

) 

419 

54 

) 

420 

85 

+ 

421 

53 

( 

422 

93 

= 

423 

00 

0 

424 

00 

0 

425 

06 

6 

426 

02 

2 

427 

00 

0 

428 

05 

5 

429 

55 

f 

430 

53 

( 

431 

43 

RCL 

432 

12 

12 

433 

65 

X 

434 

53 

( 
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435 

43 

RCL 

436 

11 

11 

437 

33 

x2 

438 

54 

) 

439 

54 

) 

440 

54 

) 

441 

54 

) 

442 

42 

STO 

443 

08 

08 

444 

99 

PRT 

445 

91 

R/S 

446 

76 

LBL 

447 

14 

D 

448 

53 

( 

449 

53 

( 

450 

04 

4 

451 

65 

X 

452 

89 

7T 

453 

65 

X 

454 

43 

RCL 

455 

11 

11 

456 

54 

) 

457 

65 

X 

458 

53 

( 

459 

53 

( 

460 

43 

RCL 

461 

08 

08 

462 

65 

X 

463 

01 

1 

464 

05 

5 

465 

00 

0 

466 

00 

0 

467 

00 

0 

468 

00 

0 

469 

55 

T 

470 

53 

( 

471 

89 

7T 

472 

65 

X 

473 

02 

2 

474 

65 

X 

475 

43 

RCL 

476 

12 

12 

477 

54 

) 

478 

54 

) 

479 

54 

) 

480 

55 

T 

481 

53 

( 

482 

53 

( 

483 

53 

( 

484 

53 

( 

485 

53 

( 

486 

43 

RCL 

487 

07 

07 

488 

55 

T 

489 

43 

RCL 

490 

12 

12 

491 

54 

) 

492 

33 

xz 

493 

54 

) 

494 

75 

- 

495 

01 

1 

496 

54 

) 

497 

33 

X2 

498 

54 

) 

499 

85 

+ 

500 

53 

( 

501 

43 

RCL 

502 

08 

08 

503 

33 

x2 

504 

54 

) 

505 

54 

) 

506 

54 

) 

507 

42 

STO 

508 

09 

09 

509 

99 

PRT 

510 

91 

R/S 

511 

76 

LBL 

512 

15 

E 

513 

53 

( 

514 

53 

( 

515 

04 

4 

516 

65 

X 

517 

89 

ir 

518 

65 

X 

519 

53 

( 

520 

43 

RCL 

521 

11 

11 

522 

33 

x2 

523 

54 

) 

524 

54 

) 

525 

55 

T 

526 

53 

( 

527 

53 

( 

528 

53 

( 

529 

53 

( 

530 

53 

( 

531 

43 

RCL 

532 

07 

07 

533 

55 

T 

534 

43 

RCL 

535 

12 

12 

536 

54 

) 

537 

33 

x2 

538 

54 

) 

539 

75 

- 

540 

01 

1 

541 

54 

) 

542 

33 

X2 

543 

53 

) 

544 

85 

+ 

545 

53 

( 

546 

43 

RCL 

547 

08 

08 

548 

33 

X2 

549 

54 

) 

550 

54 

) 

551 

54 

) 

552 

42 

STO 

553 

10 

10 

554 

99 

PRT 

555 

66 

PAU 

556 

53 

( 

557 

43 

RCL 

558 

09 

09 

559 

75 

- 

560 

43 

RCL 

561 

10 

10 

562 

54 

) 

563 

42 

STO 

564 

11 

11 

565 

99 

PRT 

566 

91 

R/S 

120 


APPENDIX  D 
NUMERICAL  INTEGRATION  TI  59  PROGRAM 

A.   INTRODUCTION 

In  order  to  perform  the  numerical  integration  of  an (a) da, 
a  standard  Texas  Instrument's  program  was  used  [Ref.  25]. 
This  program  performs  the  integration  by  using  Simpson's 
discrete  approximation  based  on  the  following  expression 

x  . 

I  =  /  n  f(x)dx  «  ;j(f  +4f !+2f2+4f 3+2f4  + +2fn— +4fn  x  +  fn) 

xo 


where 


f(x)  must  be  known  at  n+1  equally  spaced  points  (f  -f  ) . 

x  -x 
*  =  "V*  »  *n>Xo 

n+7  -  number  of  data  registers  available 
n = number  of  subintervals  =  2,  4,  6,  .... 


B.   PROGRAM  INSTRUCTIONS 

The  program  is  taken  from  the  master  library  program 
package  by  using  the  code  2nd  Pgm  10  on  the  calculator. 
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The  user's  instructions  are  as  follows 


Procedure 

Enter 

Press 
A 

Display 

Enter  data 

Subintervals 

1 
n 

Enter  data 

h 

B 

h1 

Enter  data 

Function  values 

0 

C 

0 

fo 

R/S 

^ 

fl 

R/S 

h1 

f  R/S        f  1 

n  '         n 

Calculate  -  D         I1 


These  values  are  printed  automatically  if  the  calculator 
is  connected  to  PC-100A  Print  Cradle. 


The  program  steps  are  listed  below  giving  location  (LOC) , 
key  symbol  (KEY) ,  and  comments. 


000 

76 

LBL 

001 

11 

A 

002 

53 

( 

003 

50 

lxl 

004 

42 

STO 

005 

05 

05 

006 

55 

T 

007 

02 

2 

008 

54 

) 

009 

42 

STO 

010 

02 

02 

011 

22 

INV 

012 

59 

INT 

013 

29 

CP 

014 

22 

INV 

015 

67 

EQ 

016 

52 

EE 

017 

43 

RCL 

018 

05 

05 

019 

99 

PRT 

020 

92 

RTN 

021 

76 

LBL 

022 

52 

EE 

023 

00 

0 

024 

35 

1/x 

025 

92 

RTN 

026 

76 

LBL 

027 

12 

B 

028 

42 

STO 

029 

03 

03 

030 

99 

PRT 

031 

92 

RTN 

032 

76 

LBL 

033 

13 

C 

034 

53 

( 

035 

24 

CE 

036 

85 

+ 

037 

32 

x  <  t 

038 

06 

6 

039 

54 

) 

040 

42 

STO 

041 

01 

01 

042 

32 

x^t 

043 

98 

ADV 

044 

92 

RTN 

045 

76 

LBL 

046 

50 

lxl 

047 

76 

ST* 

048 

01 

01 

049 

32 

x<t 

050 

01 

1 

051 

44 

SUM 

052 

01 

01 

053 

32 

x£t 

054 

99 

PRT 

055 

92 

RTN 

056 

61 

GTO 

122 


057 

50  lxl 

058 

76  LBL 

059 

14  D 

060 

53  ( 

061 

43  RCL 

062 

05  05 

063 

85  + 

064 

06  6 

065 

54  ) 

066 

42  STO 

067 

01  01 

068 

7  3  RC* 

069 

01  01 

070 

42  STO 

071 

04  04 

072 

7  6  LBL 

073 

45  yx 

074 

01  1 

075 

22  INV 

076 

44  SUM 

077 

01  01 

078 

53  ( 

079 

7  3  RC* 

080 

01  01 

081 

65  x 

082 

04  4 

083 

54  ) 

084 

44  SUM 

085 

04  04 

086 

01  1 

087 

22  INV 

088 

44  SUM 

089 

01  01 

090 

22  INV 

091 

97  DSZ 

092 

02  02 

093 

33  x2 

094 

53  ( 

095 

7  3  RC* 

096 

01  01 

097 

65  x 

098 

02  2 

099 

54  ) 

100 

44  SUM 

101 

04  04 

102 

61  GTO 

103 

45  yx 

104 

76  LBL 

105 

33  x2 

106 

7  3  RC* 

107 

01  01 

108 

44 

SUM 

109 

04 

04 

110 

53 

( 

111 

43 

RCL 

112 

03 

03 

113 

55 

T 

114 

03 

3 

115 

54 

) 

116 

49 

PRD 

117 

04 

04 

118 

43 

RCL 

119 

04 

04 

120 

98 

ADV 

121 

99 

PRT 

122 

92 

RTN 
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APPENDIX  E 
TURN  RATE  LIMITATION  TI  59  PROGRAMS 

A.  INTRODUCTION 

The  turn  rate  limitation  calculations  are  divided  into 
two  parts  with  separate  programs. 

B.  PART  I  PROGRAM 

For  a  given  initial  detection  range  and  a  given  ship 
speed,  Part  I  performs  a  sequential  calculation  of  cj>  ,K  , 
$_.  ,K_, ,  $*  and  K*  for  each  of  the  maximum  turn  rates  investi- 
gated,  together  with  the  "no  limit"  conditions  (V  )    and 


s'NL 


K^T    based  on 


.    -1    ,  3pl<i>|max      .Qr>o 

$A  =   sin        (      v •  ^90 

s 


*B   =    180°    -    0A 

V(l+cos   (J>  R)P 

A'B  '      | 6|  (sin  *A    B)P"2 

1     'max  XA,B 


<j>*   =   cos_1(P-2) 


V    (1+P/2)P[1-(P/2)2]1"P/2 
1 Y ' max 
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For  p>2,  the  solution  for   *  and  K*  are  not  valid. 
Furthermore,  for  Vg  (Vs)NL  no  solutions  are  valid  for  any 
of  the  quantities. 

The  program  is  based  on  the  fixed  torpedo  speed  of 
VT  =  35  kts  (18  m/s)  and  an  initial  detection  range  of 
Tq  =  =  R,,  =  250  m. 

The  user  instructions  for  the  program  are  as  follows: 


Prodcedure 


Enter 


Press 


Display 


Enter  data 

Enter  data 

Calculate  <$>* 
and  K* 

Calculate  d>„ 


and  0 


3 


Calculate  K, 
and 


*b 


Maximum  Turn  Rate 
Ship  speed 


A 
B 


D 
E 


I* 


max 


*  V  *" 


<f>*,K 


A'  B 


KA'KB 


"These  values  are  printed  automatically  if  the  calculator 
is  connected  to  PC-100A  Print  Cradle. 


A  listing  of  the  program  steps  follows  giving  location 
(LOC) ,  code  (COD) ,  key  symbol  (KEY) ,  and  comments. 


000 

76 

L3L 

001 

11 

A 

002 

42 

STO 

003 

01 

01 

004 

99 

PRT 

005 

91 

R/S 

006 

76 

LBL 

007 

12 

B 

008 

42 

STO 

009 

02 

02 

010 

99 

PRT 

011 

25 

CLR 

012 

69 

OP 

013 

00 

00 

014 

03 

3 

015 

03 

3 

016 

69 

OP 

017 

04 

04 

018 

53 

( 

019 

01 

1 

020 

08 

8 

021 

55 

T 

022 

43 

RCL 

023 

02 

02 

024 

54 

) 

025 

42 

STO 

026 

03 

03 

027 

69 

OP 

028 

06 

06 

029 

25 

CLR 

030 

69 

OP 

031 

00 

00 

032 

02 

2 

033 

06 

6 

034 

03 

3 

035 

01 

1 

125 


036 

02 

2 

037 

07 

7 

038 

69 

OP 

039 

04 

04 

040 

53 

( 

041 

43 

RCL 

042 

03 

03 

043 

65 

X 

044 

03 

3 

045 

54 

) 

046 

42 

STO 

047 

10 

10 

048 

69 

OP 

049 

06 

06 

050 

25 

CLR 

051 

69 

OP 

052 

00 

00 

053 

04 

4 

054 

02 

2 

055 

03 

3 

056 

06 

6 

057 

03 

3 

058 

01 

1 

059 

02 

2 

060 

07 

7 

061 

69 

OP 

062 

04 

04 

063 

53 

( 

064 

53 

( 

065 

53 

( 

066 

03 

3 

067 

65 

X 

068 

01 

1 

069 

08 

8 

070 

54 

) 

071 

65 

X 

072 

43 

RCL 

073 

01 

01 

074 

54 

) 

075 

34 

/x 

076 

54 

) 

077 

42 

STO 

078 

11 

11 

079 

69 

OP 

080 

06 

06 

081 

91 

R/S 

082 

76 

L3L 

083 

13 

C 

084 

25 

CLR 

085 

69 

OP 

086 

00 

00 

087 

03 

3 

088 

03 

3 

089 

02 

2 

090 

03 

3 

091 

05 

5 

092 

01 

1 

093 

69 

OP 

094 

04 

04 

095 

53 

( 

096 

53 

( 

097 

53 

( 

098 

43 

RCL 

099 

03 

03 

100 

55 

T 

101 

02 

2 

102 

54 

) 

103 

22 

INV 

104 

39 

COS 

105 

54 

) 

106 

54 

) 

107 

42 

STO 

108 

04 

04 

109 

69 

OP 

110 

06 

06 

111 

25 

CLR 

112 

69 

OP 

113 

00 

00 

114 

02 

2 

115 

06 

6 

116 

05 

5 

117 

01 

1 

118 

69 

OP 

119 

04 

04 

120 

53 

( 

121 

53 

( 

122 

53 

( 

123 

53 

( 

124 

43 

RCL 

125 

02 

02 

126 

65 

X 

127 

53 

( 

128 

53 

( 

129 

01 

1 

130 

85 

+ 

131 

53 

( 

132 

43 

RCL 

133 

03 

03 

134 

55 

T 

135 

02 

2 

136 

54 

) 

137 

54 

) 

138 

45 

X 

y 

139 

43 

RCL 

140 

03 

03 

141 

54 

) 

142 

54 

) 

143 

65 

X 

144 

53 

( 

145 

53 

( 

146 

01 

1 

147 

75 

- 

148 

53 

( 

149 

53 

( 

150 

43 

RCL 

151 

03 

03 

152 

55 

T 

153 

02 

2 

154 

54 
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C.   PART  II  PROGRAM 

Given  K  from  Part  I,  Part  II  program  performs  the  probe 
calculations  of  the  corresponding  initial  approach  angles 
^o*A'  ^o^B'  ^o^  *'  and  ^o^NL  based  on  the  relationship 


f(*o»  "  Ka.b^nl^^o'^b,*,^"1 

-ro[(l+(cos*o)AjB^NL)P]  =  0 


The  program  used  a  fixed  torpedo  speed  of  V_  =  35  kts  (18  m/s) 
and  an  initial  detection  range  of  r  =  R^  =  250  m.   The  purpose 
of  this  general  probe  program  is  to  locate  roots  of  the  given 
function  y  =  f  (cj>  )  to  evaluate  the  slope  of  the  tangent  line, 
and  to  find  the  maximum  and  minimum  points  on  a  graph.   We 
will  only  use  the  first  feature. 
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The  program  requires  a  subroutine  for  the  function  to 
be  investigated.   This  subroutine  starts  at  program  location 
140  and  is  located  at  label  2nd  C.   The  user's  instructions 
for  the  program  are  as  follows: 


Procedure 


Enter 


Press 


Enter  data 

Starting  Value 

of 

*o 

2nd 

D 

*o 

Enter  data 

The 

increment  Ad> 
o 

2nd 

E 

A*o 

Calculate 

f(VA*Q) 

— 

A 

f(VA*o) 

Calculate 

f(*o-a*0) 

— 

B 

f(*0-A*Q) 

Display  current 

value  of  $ 
o 

E 

*o 

If  the  value  of  A<J>   is  chosen  too  large,  A<(>   may  be 
replaced  by  Ac?  /10  by  pressing  label  2nd  A.   Similarly,  if 
a  larger  value  of  A<J>   is  required,  A<J)   can  be  replaced  by 
10A<{)   by  pressing  label  2nd  B.   A  listing  of  the  program 
steps  follows,  giving  location  (LOC) ,  code  (COD),  key  symbol 
(KEY),  and  comments. 

Associated  with  the  subroutine,  it  should  be  noted  that 
the  value  of  K,r   and  p  are  entered  separately  in  the 
memory  locations  10,  11,  and  12,  respectively. 
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Fig.  2.   Location  of  Weather  Stations  Along  the  Norwegian 
Coast. 
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Fig.  4.   Histogram  of  Significant  Wave  in  Percent  per  Year 
from  the  Weather  Station  "Andenes." 
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Fig.  5.   Resonant  Bubble  Density  at  12  kHz  as  a  Function  of 
Depth. 

136 


10000 
8 

6 


'  Resonansfrekvens  38  kHz 

Resonance   frefjuency  t38  kHz 


•    9.0  knt. 
a  12.0  - 


0    1     2    3    U    5    6    7    8    9    10  11    12  13  14  15  16 
Depth  (r.)  Qybde  i  [ml 


Fig.  6.   Resonant  Bubble  Density  at  38  kHz  as  a  Function  of 
Depth. 


137 


10 


8   Resonans  frekvens    120kHz 

£  |Resp.nai^ce  .fceq^jency  IgO  kHz.^_ 


8 

So 


B 

.0 


2  - 

6 

10  — 


2 — 


■    5 

8 
6 


l 


c 

0) 

a 


a 
a 


3     4 


PQ 


4~ 


2  — 


10 
8 
6 

4— 

fr- 
it? 


•  12  knt. 
X35  » 
a  40" 
a  45" 


0    12    3 


4    5    6    7    8    9  10  11   12  13  14  15 
Depth  (n)  Dybde  i  [ml 


Fig.  7.   Resonant  Bubble  Density  at  120  kHz  as  a  Function 
of  Depth. 

138 


LLmtti    of    prevailing  noise 
■  Wind- dependent  bubble  dnd  s pray  noise 

Heavy    precLpitaUon 

Heavy    traffi 


Frequency  (kHz) 


Fig.  8.   Ambient  Noise  Level  Curves. 
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Fig.  9.   Typical  Sound  Speed  Profiles  in  Norwegian  Coastal 
Waters. 
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Fig.  12.   Frequencies  of  Wind  Directions  in  Percent  for 

Stations  Along  the  Coast  from  "Hillesoy"  to  "Ona." 
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Fig.  13.   Frequencies  of  Wind  Directions  in  Percent  for 

Stations  Along  the  Coast  from  "Myken"  to  "Furuholmen. " 
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Fig.    14.      Absorption  Coefficient   in  dB/m  as   a  Function  of 
Frequency. 
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Fig.  15.   Average  Radiated  Spectrum  Level  for  Surface  Ship 
as  a  Function  of  Speed  in  kts. 
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Fig.  17.   Spectrum  Level  as  a  Function  of  Total  Gas  Content 
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Fig.  20.   Required  Input  S/N  Ratio  and  BT  Product  for  Various 
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Fig.  21.   Probability  of  Detection  Versus  Input  S/N  Ratio  for 
Various  Thresholds. 
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Fig.    23.       - 


(gp)    ssoi   buxpeaads    XHOT^aqds 
201ogR  -aR  Versus   R  for   a  Frequency  of   30   kHz. 
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Fig.    24a.      Geometry   at   the   Sea   Surface   Scattering, 
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Fig.  24b.   Specular  Scattering  Geometry. 
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Fig.  25.   Rough  Surface  Shadowing  Geometry 
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Given   :         h4     »  Source    depth 

hr  =  receiver  depth 
£><$>    =  haU  beam  width  of  receiver 
W  =  wind,  speed 
Varying:     R0      detection   range 
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Fig.    26.      Calculation  Scheme  of  the   Surface   Scattering  Effect, 
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Fia.  28.   Geometry  of  the  Subsurface  Propagation  Model 
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Fig.  29.   a)   Chains  of  Successive  Scattering. 

b)   Scattering  Pattern  Going  Through  the  Same 
Scatterer  More  than  Once. 
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Fig.    30.      Scattering  Pattern. 
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Fig.  34.   a  (a) n (a) da  for  60  kHz 
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Fig.  38.   Pursuit  Homing  Geometry, 
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Fig.  40.   Hit  Criterion  Geometry, 


172 


-80' 


-90' 


-100 


-10° 


15—  - 
20- 
25—  — 

l 

180° 


120" 


/ 


s 


l"30 


140 


\  160 

170° 


150' 


Fig.  41.   AOB-limitation  Versus  Ship  Speed  with  Parameter 
Maximum  Turn  Rate. 
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